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ABSTRACT 

This Final Report is a summary of the work that was performed under Contract NAS8- 
38075 between NASA Marshall Space Flight Center and Remote Sensing Systems from Sep- 
tember 1989 to September 1992. The primary accomplishment was the delivery of SSM/I 
data tapes containing sensor and geophysical products. In all, 515 tapes (80 gigabytes) were 
delivered. These tapes contained the F08 SSM/I data for the period from July 1988 
through December 1991 and the F10 SSM/I data for the period from December 1990 
through December 1991. For the F08 SSM/I, a data inventory was compiled and an en- 
gineering assessment was done. Ephemeris tables for the F08 and F10 spacecrafts were 
computed. Scientific studies on the oceanic wind vector and water vapor field were pub- 
lished, and color atlases of monthly ocean products were produced. This investigation was 
part of NASA’s WETNET program. 
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1. INTRODUCTION 

This Final Report provides a summary of the work that was performed under Contract 
NAS8- 38075 between NASA and Remote Sensing Systems (RSS). The period of perfor- 
mance for this contract was September 1989 to September 1992. A detailed description of 
the investigation can be found in a series of reports and scientific papers that were 
delivered during the course of the contract. These documents are listed below on page 9 
and are cross-referenced throughout this summary. The principal document (Document 1) 
is included herein as Appendix A. 

The main thrust of this effort was the delivery of quality controlled SSM/I data 
products, both Level- 1 sensor data and Level-2 geophysical products. This required an ex- 
tensive data processing effort in conjunction with a thorough engineering assessment of the 
SSM/I in-orbit operation. The production of accurate geophysical products required scien- 
tific investigations of the oceanic wind field and columnar water vapor. One particularly 
interesting finding was that the microwave emission from the ocean depends on the wind 
direction relative to the observation direction. With respect to water vapor, comparisons 
with radiosondes and ECMWF fields indicated that the SSM/I water vapor is an extremely 
accurate parameter that may be of considerable use to climate change studies. 

During this investigation, we have had the good fortune to interact with a large number 
of colleagues in associated disciplines via the WETNET workshops, E-mail, and the many 
phone calls about SSM/I. We consider this interaction to be an important element of this 
investigation. 


2. DELIVERY OF SSM/I DATA TAPES 

During the course of this contract, RSS delivered 336 F08 SSM/I antenna temperature 
(T A ) tapes and 94 F10 T A tapes to Marshall Space Flight Center (MSFC). A duplicate set 
of T a tapes was also delivered to the National Snow and Ice Data Center (NSIDC). The 
F08 tapes covered the period from July 1988 through December 1991, and the F10 tapes 
covered the period from December 1990 through December 1991. Approximately 6000 
DEF tapes from NOAA/NESDIS were ingested in order to produce the 430 T A tapes. 

In addition to the T A tapes, RSS delivered 84 F08 SSM/I geophysical tapes to MSFC. A 
duplicate set of 84 geophysical tapes was also delivered to the NASA Ocean Data System 
(NODS) at the Jet Propulsion Laboratory (JPL). These tapes covered the period from July 
1988 through December 1991. Software for reading the tapes was also provided along with 
comprehensive User Manuals (Documents 1 and 2). 


3. INVENTORY OF THE F08 DATA AND ERRONEOUS DATA PERIODS 

We compiled and delivered an inventory of all F08 SSM/I data (Document 4). This in- 
ventory -is a collection of figures that show time periods of missing data and erroneous data. 
There is one figure for each semimonthly period from July 1987 through December 1991. 
The figures show the differences between a 90-second average of the individual T A ’s minus 
an average 1° latitude by 1° longitude T A map for the previous semimonthly period. There 
are 16 horizontal lines in each figure, and each line corresponds to a Greenwich day. The 
90-second T A differences are shown as points along the horizontal lines. 
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Under normal conditions, the points lie near the horizontal lines, with an occasional 
small outlier. Storms, coast lines, and varying sea ice can cause larger differences, but only 
for a short period of time. The periods of erroneous data are shown as clusters of widely 
varying points. In addition to revealing periods of erroneous data, the figures provide a 
complete inventory of the SSM/I data set. They can be used to find periods of missing 
data. These data gaps are shown by the absence of plotted points. 

We visually inspected all of the figures in the F08 inventory and identified the begin 
and end times for the periods of erroneous data. A small ASCII data file called BADLOC 
containing these times was delivered on a 3.5" inch diskette. The time windows for er- 
roneous data are somewhat larger than the clusters of data outliers shown in the figures in 
order to ensure that all the erroneous data are excluded. The problem of erroneous data 
was particularly bad during the second half of 1988. Then, in January 1989 and thereafter 
the frequency of erroneous data greatly reduced. After 1988, there are, on the average, 
about one or two erroneous data periods each month. 

There appears to be two types of erroneous data. First, the observation time is some- 
times in error, and as a result the geolocation routine mislocates the observations. In this 
case the T A values are correct but the latitudes and longitudes are in error. The second 
type of erroneous data is the occurrence of corrupt T A values that are obviously incorrect. 
The cause of the corrupt T A ’s is still not known. Both types of erroneous data usually oc- 
cur over a complete DEF data file and hence are detected using the above procedure. 
When doing statistical analyses of the SSM/I data, it is very important to exclude these time 
periods of erroneous data. Otherwise, the statistics can be significantly distorted. 


4. ENGINEERING ASSESSMENT OF SSM/I OPERATION 

The in-orbit operation of the F08 SSM/I from July 1987 through December 1991 was 
analyzed in detail. Time series of the following radiometer parameters were produced: 

1. Thermistor readings of the warm load and radiometer receivers. 

2. The rms variation of the SSM/I cold-space observations 

3. The rms variation of the SSM/I warm-load observations 

4. The radiometer noise temperature 

5. The radiometer gain 

The time series are shown by a collection of color-coded figures, with different colors 
showing the 7 SSM/I channels. In addition to these figures, logs were maintained of the 
automatic gain control settings and the number of anomalous T A observations. 

The time series show the stability (and in some cases the instability) of the SSM/I chan- 
nels. For example, the figure of thermistor readings shows a very large heating cycle with 
a peak between 30K and 60K. above normal, depending on component, every winter due to 
increased solar illumination on the SSM/I. 

During the first heating cycle in December 1987, the 85V temperature resolution begins 
to degrade. After the heating cycle is over, the 85V channel continues to degrade. After 
April 1988, the noise in the 85V channel exceeds 20 K. Between May 1988 and January 
1989, there are time windows in which the noise subsides. The possibility exists that useful 
information can be obtained from the 85V channel during these time windows. After 
January 1989, the noise becomes very large, and the channel is not usable. 
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The 85H channel shows a slight degradation in the temperature resolution during the 
first heating cycle, but then recovers. The second heating cycle in December 1988 seems to 
cause a small but permanent degradation in the 85H channel. Then, the third heating cycle 
in December 1989 seriously damages the 85H channel, and the noise increases to between 5 
and 10 K. Curiously, during the fourth heating cycle in December 1990, the 85H channel 
temporarily recovers, but after January 1991 the noise becomes very large, exceeding 20 K. 

The other channels also experience a small increase in noise during the heating cycles. 
Fortunately, the lower frequency channels are not permanently affected by the heating, and 
noise level returns to normal after the heating cycle is over. We also verified that the 
radiometer noise computed from the cold and warm counts is in good agreement with that 
computed from the radiometer noise temperature and bandwidth. 

In addition to analyzing the time series for radiometer noise, we also investigated the 
variation of the radiometer gain with time. Under normal operating conditions, the varia- 
tion in the radiometer gain G r is primarily due to changes in the temperature of the 
receiver. The large heating cycle experienced by the F08 SSM/I provided a means to in- 
vestigate the variation of G r with temperature. The gain of the 19V and 19H channels 
showed little dependence on temperature, while the gain of the 22V channel increased with 
increasing temperature. In contrast, the gain of the 37V and 37H channels decreased with 
increasing temperature. The 85V and 85H gains also decreased with temperature, and then 
tended to zero during the degradation phase. 

The dependence of G r on temperature is one of the primary reasons for doing an exter- 
nal calibration each scan. The temperature variation experienced during a single orbit will 
produce changes in G r . Unless frequent calibrations are performed, the changing gain will 
introduce errors in the computation of the antenna temperatures. We verified that the 
derivative of G r with respect to temperature found during the winter heating cycle (i.e., 
long-time scale) is the same as that obtained from an analysis of the gain versus tempera- 
ture variation experienced within individual orbits (short-time scale). 

We monitored changes in the automatic gain control (AGC) setting during the 1987-1991 
period. The purpose of the AGC is to compensate for the G r variation so that the counts 
stay within the range of the digital converter. There are 16 AGC settings, and each setting 
adds 0.5 dB attenuation. The AGC system seemed to function properly. 

As part of the engineering assessment, we also found a systematic error in the SSM/I 
T a ’s near one edge of the swath. The T A ’s show a roll-off of about 1 K near the scan 
edge. One explanation for this is that the feedhorn begins to partially see the cold-sky 
reflector near the end of the Earth-viewing scan. A correction table was produced in or- 
der to remove this systematic error. 

The engineering assessment is further described in Document 1. 
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5. GEOLOCATION PROCESSING AND EPHEMERIS TABLES 


While processing the 1988 SSM/I data, we found a number of problems with the 
spacecraft and Earth locations that were being computed at Fleet Numeric Oceanography 
Center (FNOC). The cell latitudes and longitudes on the DEF tapes produced by FNOC 
were sometimes in error due to the following problems: 

1. There are some algorithm errors in the FNOC data processing software. 

2. The satellite ephemeris is sometimes incorrect due to spacecraft tracking errors and orbit 
prediction errors. This problem is particularly severe during periods of increased solar ac- 
tivity. 

3. The boresight nadir angle and the alignment of the SSM/I instrument relative to the 
spacecraft are slightly misspecified. 

There have been a number of FNOC algorithm errors. Sometimes the wrong Julian day 
is assigned to the orbit. Occasionally, the time for the spacecraft location and the scan time 
are very different. The two times should be within 60 s of each other, but at the begin 
and end of an FNOC DEF file, the two times can be as much as 40,000 s apart. The 
spacecraft longitude is sometimes out-of-bounds. There were also other problems that have 
now been corrected at FNOC. For example, the ephemeris was not updated frequently 
enough, and the location of cell 128 was in error. 

In response to these problems, we developed our own routine for computing the latitudes 
and longitudes rather than using the DEF values. The input to our geolocation routine is 
the satellite ephemeris for a ±7 day period centered on the orbit being processed. The 
ephemeris is first subjected to quality control to remove any data that are clearly erroneous. 
Next the ephemeris is smoothed to remove the noise. The smoothed ephemeris is then used 
to compute the SSM/I cell latitudes and longitudes. 


j.his procedure of computing latitude and longitude from a smoothed ephemeris was in- 
itiated with the January 1, 1989 F08 T A tape. Prior to 1989, the latitudes and longitudes on 
the T A tapes are the same as the values on the DEF tapes. Fortunately, an analysis of the 
ephemeris prior to 1989 does not show any large errors except for a brief period during the 
second half of February 1988. For these pre-1989 tapes, subroutine DECODE does the 
sensor alignment and the cell 128 correction while reading the tape. 

Our geolocation routine is free of the algorithm errors listed above except for the Julian 
day error. As each orbit is processed, we perform a number of cross-checks to verify the 
Julian day. Usually these cross-checks detect the Julian day error, which can then be cor- 
rected. However, a few of the Julian day errors will go undetected, and these ‘erroneous’ 
orbits will be written to the T A tape. These erroneous orbits are later detected by the 
quality control routine discussed in Section 3 above. 

The third geolocation problem listed above is the misalignment of the SSM/I sensor rela- 
tive to the spacecraft and the misspecification of the boresight nadir angle. For F10, this 
problem is corrected by our geolocation routine, and the locations on the T A tapes are cor- 
rect. For F08, part of the correction is done by our geolocation routine and part is done by 
DECODE, as is explained in Document 1. 
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An ephemeris table for the F08 spacecraft was produced. This table gives the value of 
six orbit parameters for every F08 orbit from July 1987 through December 1991. The six 
orbit parameters are the following: 

1. orbit period (seconds) 

2. local time of the ascending node (seconds) 

3. inclination angle (degrees) 

4. major semi-axis (km) 

5. eccentricity 

6. angle of perigee (degrees) 

Figures were produced that show the variation of these six parameters with time. The 
period P and the major semi-axis a show a slow, steady decrease with time. For example, 
at the beginning of the F08 mission in July 1987 the values for P and a were 6118.5 s and 
7229.3 km. Four years later, the value are 6111.0 s and 7223.3 km, respectively. The as- 
cending node local time t L steadily increases with time. For F08 the change is small, being 
about 70 s/year. The inclination angle a is slowly decreasing with time at a rate of about 
-O.Or/year for F08. The eccentricity e and angle of perigee T 0 are periodic. The period 
is about 1800 orbits and the amplitude remains constant over time. Note that the period 
and the major semi-axis are redundant parameters in the sense that they are related by 
Kepler’s third law of motion. The expression 21.6106 P 273 yields values for a that are 
within 0.2 km of the values in the orbit parameter table. 


6. EVALUATION OF NOAA’S LEVEL IB SSM/I PRODUCT 

We obtained NOAA’s new Level IB SSM/I tape product. Unfortunately, there was very 
little documentation to support the Level IB product. We compared the Level- IB SSM/I 
data with the data on the FNOC DEF tapes. The scan times and cell latitude and lon- 
gitudes were different. Due to the lack of documentation, we were unable to resolve the 
differences. Also, the Level- IB product contained no information on the satellite location 
time, and the satellite location vectors were not defined in the documentation. For the 
Level- IB product, some type of averaging had been done to the hot load thermistor counts. 
Again, the documentation did not provide any information on how the averaging was done. 


7. ANALYSIS OF SSM/I WIND SPEED AND WIND DIRECTION SIGNAL 

We investigated the possibility of retrieving both wind speed and direction from SSM/I. 
The SSM/I observations were collocated with buoy reports from the National Data Buoy 
Center (NDBC). A physically- based algorithm was used to retrieve the wind speed. The 
rms difference between the SSM/I and buoy wind speed was 1.6 m/s for 3321 comparisons. 
We found that the SSM/I minus buoy wind speed difference is correlated with wind direc- 
tion. This correlation is due to a wind direction signal in the brightness temperatures (T B ). 
For moderate wind speeds, the v-pol upwind T B is about 2 K higher than the downwind 
T b - I"° r h-pol observations, the T B maximum occurs at crosswind rather than upwind and 
is approximately 3.5 K higher than the minimum at downwind. When this wind direction 
signal is removed, the rms difference between the SSM/I and buoy winds reduces to 1.3 
m/s. The wind direction signal was used to make global, low- resolution maps of the 
monthly mean oceanic wind vector. Finally, we assessed the wind direction sensing 
capability of a prospective two-look satellite radiometer. Preliminary computer simulations 
indicate a direction accuracy of 15°. The results of this investigation have been published 
in IEEE Geoscience and Remote Sensing (Document 5). 
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8. ANALYSIS OF SSM/I WATER VAPOR PRODUCT 

The SSM/I water vapors were evaluated by comparing them to radiosonde observations 
and to ECMWF fields. The rms difference between the SSM/I columnar water vapor and 
the radiosonde value is 3 mm for hourly comparisons and 4 mm for monthly comparisons. 
The SSM/I vapor fields agree with the ECMWF fields over most ocean areas, but there are 
significant differences in the dry air masses over the eastern tropical oceans, particularly in 
the southern hemisphere. Radiosonde reports in these areas indicate that the ECMWF 
vapor fields are too high, and the SSM/I values seem accurate. The results of this inves- 
tigation have been published in JGR (Document 6). 


9. ATLAS OF MONTHLY OCEAN PRODUCTS 

We produced a 4.5 year data set of ocean products derived from the F08 SSM/I. The 
data set covers the period from July 1987 through December 1991. The products include 
the near-surface wind speed W, the columnar atmospheric water vapor V, and the columnar 
atmospheric liquid water L (combination of cloud and rain water). The ocean products 
have been averaged into monthly IMatitude by IMongitude bins. these monthly 
averages are put on a single 6250-bpi, 2400-ft tape, called the "Monthly Ocean Tape". This 
tape contains 53 files, which correspond to each month from July 1987 through December 
1991, except that December 1987 is missing because the SSM/I was turned off. Note that 
in July 1987 and January 1988, the SSM/I was turned off part of the time. The operation 
period for July 1987 is from July 9 through July 31. The operation period for January 
1988 is from January 13 through January 31. Also, starting on December 19, 1991, the 
Fleet Numerical Oceanography Center (FNOC) began limited processing of the F08 SSM/I 
data. Typically, only one orbit per day of F08 data is now being processed. The monthly 
ocean products are further described in Document 3. 

We collaborated with other investigators to merge the monthly SSM/I ocean products 
with the monthly mean distributions of the GEOSAT sea surface height, the AVHRR sea 
surface temperature, and the ECMWF surface wind components. Atlases that show these 
products as global colored maps were produced for 1987, 1988, and 1989 (Documents 7, 8 
and 9). 


10. GEOPHYSICAL SOFTWARE 

The software for computing geophysical parameters (wind speed, water vapor, and liquid 
water) from the SSM/I T^’s was delivered. This software is a set of Fortran-callable 
routines. The software computes wind speed W, water vapor V, and liquid water L by 
solving the following three equations in three unknowns: 


Tb22v = F 22 v(W,V,L) 

T B 37v = F 3 7v(W,V,L) 

T B37h = F 37 h(W,V,L) 

where T B22v , T B37v , an< * T B3 7 h are the SSM/I brightness temperature measurements for the 
following three channels: 22-GHz v-pol, 37-GHz v-pol, and 37-GHz h-pol. The unknowns 
are wind speed W (m/s), columnar water vapor V (g/cm 2 ), and columnar liquid water L 
(mg/cm 2 ). The brightness temperature model functions for the three channels are denoted 
by F 22v , F 37v , and F 37h . These T B functions have the form 

F(W,V,L) = T Bu + rtET s + (1 - E)(l + wW)(T Bd + rT^)] 
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where the channel subscript is now implicit. The quantities T Bu and T Bd are the upwelling 
and downwelling brightness temperatures due solely to the atmospheric emission and ab- 
sorption. The quantity T^. is the brightness temperature of cold space and equals 2.8 K at 
37 GHz. The emissivity of the sea surface is denoted by E. The sea-surface temperature 
is denoted by T s and is in units of Kelvin degrees. The term (1 + wW) accounts for a small 
increase in the reflected atmospheric radiation due to surface scattering (as opposed to 
simple reflection). The term r is the atmospheric transmission and is given below as a 
function of V and L. 

We assume that there is no significant scattering of radiation by the atmosphere. At 22 
and 37 GHz, scattering occurs only when there is appreciable rain or ice in the field of 
view. The upward and downward atmospheric brightness temperatures for a non-scattering 
atmosphere are given by 


T Bu = (1 - r)T u 
Tj3d = (1 - r)T d 

where T u and T d are the effective air temperatures. If the air temperature T a were con- 
stant (i.e., a zero lapse rate), then T u and T d would exactly equal T a . In reality the air 
temperature in the troposphere usually decreases with altitude, and T u and T d are weighted 
averages of the air temperature, where the weight is approximately equal to the absorption 
along the propagation path. 

The total absorption along the SSM/I viewing path is equal to ln(l/r) and is the sum of 
three components: oxygen, water vapor, and liquid water (i.e., both cloud and rain water). 
It is convenient to normalize the absorption in terms of a vertically integrated quantity 
rather than a viewing- path integrated quantity. In this way, the dependence on incidence 
angle d { is removed. For incidence angles below 60°, the ratio of the viewing-path length 
through the troposphere to the height of the troposphere is simply sec d- The total ab- 
sorption is then given by 


ln(l/r) = (A 0 + A v + A^) sec 6 i 

where A 0 , A v , and A L are the vertically integrated, or columnar, absorption components 
due to oxygen, water vapor, and liquid water, respectively. Absorption A v is proportional 
to V and A L is proportional to L. 

Document 5 gives values for E, w, T u , T d , A 0 A v , and A L . Once these values are 
specified, the T B model functions and their first derivatives with respect to W, V, and L 
can be computed. 

Having specified F(W,V,L), the three T B model function equations are simultaneously 
solved by Newton’s method. The function F(W,V,L) is quasi-linear in terms of W, V, and 
L, and the following approximation is used: 

F(W,V,L) » F(W 0 ,V 0 ,L 0 ) + (3F/dW)(W - W 0 ) + (dF/dV)(V - V 0 ) + (5F/3L)(L - L 0 ) 

The terms W 0 , V 0 , and L 0 are ‘first-guess’ values. Good choices are 8 m/s, 3 g/cm 2 , and 
10 mg/cm 2 , respectively, which represent typical values over the ocean. The final solution 
is independent of the first guess values, but the number of iterations is minimized by start- 
ing with typical values. The derivatives dF/3W, dF/dV, and dF/dL are found by analyti- 
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cally differentiating the T B model functions evaluated at (W 0 , V D , L 0 ). The above equation 
represents three linear equations in three unknowns: W, V, and L. The system of linear 
equations is inverted to yield the unknowns. 

A second iteration is then done for which the first guess values are replaced by the 
values coming from the first iteration. This procedure is repeated until the the values for 
W, V, and L converge. Typically, convergence is obtained after 3 to 5 iterations. If con- 
vergence is not obtained in 10 iterations, then the brightness temperatures are flagged as 
anomalous. For SSM/I, non-convergence occurs very rarely (less than 0.01% of the obser- 
vations) and indicates either erroneous brightness temperatures or possibly heavy rain. 

A full description of the exact algorithm is given in Document 5. 


11. PARTICIPATION IN WETNET, HDAAWG, AND SSM/I PATHFINDER 

As part of this contract, Frank J. Wentz actively participated in three NASA programs: 
WETNET, HDAAWG, and the SSM/I Pathfinder. He was the group leader for the WET- 
NET vapor and clouds groups. Participation in these programs required the preparation of 
presentations and travel. The following trips were made: 

1. A 3-day WETNET meeting at MSFC on 2/21/90 through 2/23/90 

2. Annual WETNET meeting in Santa Barbara on 2/19/91 through 2/22/91 

3. NOAA Level- IB meeting at GSFC on 2/24/91 through 2/26/91 

4. WETNET training session at MSFC on 5/28/91 through 5/30/91 

5. Annual WETNET meeting at Tallahassee on 4/20/92 through 4/23/92 

6. HDAAWG meeting at MSFC on 6/2/92 through 6/4/92 
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SECTION 1. INTRODUCTION 


This document is the first revision of the Special Sensor Microwave/Imager (SSM/I) 
Antenna Temperature (T A ) Tapes User’s Manual. The original User’s Manual for the 
SSM/I T a tapes was issued on March 25, 1988. This first revision provides the User with 
more information on the product and implements several improvements and extensions. 
There are now three SSM/I’s in operation. The first was launched on June 19, 1987, 
aboard the Defense Meteorological Satellite Program (DMSP) F08 satellite. The second was 
launched on December 1, 1990, on the DMSP F10 satellite. And, the third was launched on 
November 28, 1991, on the Fll satellite. The Air Force plans to maintain at least two 
SSM/I’s in operation through the end of this decade. We plan to continue providing Users 
with T a tapes for this series of satellites. 

The SSM/I T a tapes are a well-documented, user-friendly product that contains all the 
information coming from the SSM/I. Essentially, the tapes contain geolocated antenna tem- 
peratures along with the sensor calibration data. For each satellite, there are eight 6250-bpi 
T a tapes per month. Software is provided to convert the antenna temperatures to bright- 
ness temperatures (Tg). The integrity of the T A tapes has been established by wide usage 
of this product for the last several years. 

The Fortran subroutine DECODE provides the interface between the T A tapes and the 
User. The input to DECODE is the string of bytes contained in a single tape block. 
DECODE organizes the data, converts the bytes to the proper engineering units, and stores 
the results in a common area accessible to the User. An ASCII file containing the 
DECODE source code is on the enclosed 3.5" diskette. The diskette also includes the file 
BADLOC, which is a list of time periods containing erroneous F08 data. 

In addition to reading the T A tapes, DECODE performs other functions such as a cor- 
rection for F08 geolocation errors, the computation of incidence angle, and the computation 
of brightness temperatures. Some of the functions depend on which satellite is being 
processed. DECODE automatically determines the satellite number and then selects the ap- 
proximate procedure. There have been two minor changes in the T A tape format (see Sec- 
tion 10). DECODE automatically adjusts to these format changes by looking at the scan 
time. We advice Users to use DECODE rather than writing their own ingest routine. If 
the User has some unusual application, we will gladly provide assistance in modifying 
DECODE to suit the application. However, the Revision- 1 DECODE should be sufficient 
for nearly all applications. 

There are several improvements and extensions of the Revision- 1 DECODE. These in- 
clude the following: 

1. DECODE is now independent of the byte-ordering of the Fortran compiler. For ex- 
ample, Fortran compilers for many mini-computers use the convention that the most- 
significant byte in a word is stored first in memory. In contrast, Fortran compilers for 
PC’s and Vax machines use the opposite convention that the least-significant byte is stored 
first. The original version of DECODE used equivalences to relate bytes to words, and 
hence it was dependent on the byte-ordering convention. (Users with PC’s and Vax 
machines had to modified DECODE.) Revision- 1 DECODE does not use equivalences. 
Rather, it directly computes the sensor data from the individual bytes and executes properly 
using either type of compiler. 
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2. Revision- 1 DECODE is setup to process multiple satellites. DECODE automatically 
determines the satellite number and then selects the appropriate procedures for geolocation 
error correction, incidence angle computation, and T A bias removal. The satellite number is 
stored in the common area /OUTDAT/ that is accessible to the User. 

3. Revision- 1 DECODE corrects the geolocation errors for the F08 SSM/I. These correc- 
tions include an along- track adjustment, a yaw adjustment, and a correction for cell 128 at 
the edge of the scan. With these corrections, the location accuracy of the SSM/I cells is 
probably about 10 km or better. For F10 and subsequent satellites, the cell locations on the 
T a tapes have already been corrected, and no action by DECODE is required. 

4. A more exact computation of the incidence angle is done. The incidence angle accuracy 
is now 0.03° as compared to 0.1° for the original version. For the F08 satellite, the in- 
cidence angle is relatively constant, being 53.1° ± 0.25°. The ± 0.25° variation is due to 
the slight eccentricity (e « 0.0015) of the F08 orbit and the oblateness of the Earth. The 
eccentricity for the F10 orbit is greater (e » 0.0085), and as a consequence the incidence 
angle variation is greater. For F10, the incidence angle ranges from 52.5° to 54.0°. This 
± 0.75° variation is large enough to introduce an error is some geophysical parameters if it 
is not taken into account. 

5. The F08 antenna temperatures are corrected for an along-scan error that occurs near the 
edge of the scan where the feedhorn partially sees the cold-sky reflector. There is also an 
option for adding an offset to the antenna temperatures in order to bring them into agree- 
ment with the Wentz [1991] radiative transfer model for the ocean. The T A corrections for 
the F10 SSM/I are currently set to zero. We will provide the Users with the correction 
coefficients as soon as the T A error analysis is completed. 

In addition to these DECODE updates, a more complete description of the SSM/I T A 
data set is given herein. We have included new sections on geolocation, erroneous data 
periods, and radiometer parameters. The table BADLOC that lists periods of erroneous F08 
data is provided on the enclosed 3.5" diskette (see Section 7). We also include Figure 9 in 
Section 8 that shows when the F08 85 GHz channels are usable. For more information on 
SSM/I, we refer the User to the SSM/I report compiled by Hollinger [1989]. 

For those Users who are not interested in the details of the SSM/I data processing and 
sensor performance, we suggest that you at least read Sections 7, 9, and 11. It is important 
to excluded periods of erroneous data (Section 7), to know the organization of the T A tapes 
(Section 9), and to understand the execution of the subroutine DECODE (Section 1 1 ). 
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SECTION 2. SENSOR DESCRIPTION 


The orbits for F08 and F10 are near-circular, sun-synchronous, and near-polar, with 
an inclination of 98.8°. The altitude for F08 is 860 ± 25 km, and it is 805 ± 72 km for 
F10. The variation in altitude is due to the eccentricity of the orbit and the oblateness of 
the Earth. The F08 and F10 orbital periods are 102 and 101 minutes, respectively. The 
local time for the ascending equatorial crossing for F08 is 6:12 am. Due to a malfunction 
during launch, F10 is not exactly sun-synchronous, and its local equatorial crossing time is 
increasing at a rate of 47 minutes/year. On January 1, 1991, the F10 local equatorial cross- 
ing time was 7:39 pm. More information on the orbit parameters is given in Section 6. 

The SSM/I consists of 7 separate total-power radiometers, each simultaneously measur- 
ing the microwave emission coming from the Earth and the intervening atmosphere. Table 
1 gives the frequencies, polarizations, and temporal and spatial resolution of the 7 channels. 
Dual-polarization measurements are taken at 19.35, 37.0, and 85.5 GHz, and only vertical 
polarization is observed at the 22.235 GHz water vapor channel. 

The SSM/I uses an offset parabolic reflector of dimensions 61 by 66 cm to collect the 
microwave radiation. The reflector focuses the radiation into a corrugated, broad-band, 
7-port feedhorn. The reflector and feedhorn spin as a unit about an axis parallel to nadir. 
The rotation period is 1.9 s. A cold-space reflector and a hot reference load are attached 
to the spin axis and do not rotate. The rotating feedhorn observes the fixed cold reflector 
and hot load once each scan. In this way, calibration observations are taken every scan. 

Earth observations are taken during a 102.4° segment of the rotation, as is shown in 
Figure 1. The 102.4° arc is centered on the spacecraft subtrack and corresponds to a 
1400 km wide swath on the Earth’s surface. Relative to the spacecraft velocity vector, the 
F08 SSM/I is looking aft, and the F10 SSM/I is looking forward. The 1400-km swath and 
the orbit inclination of 98.8° provide complete coverage of the Earth in two to three days, 
except for two small circular sectors of 2.4° centered on the North and South poles. 

During each scan, the 85 GHz channels are sampled i28 times over the 102.4° arc. 
The integration period for a single sample is 3.89 ms. This sampling scheme results in 128 
v-pol footprints and 128 h-pol footprints, each having an effective 3-dB spatial resolution 
of 15 km. Observations at the lower three frequencies are only taken every other scan. 
Scans during which the lower channels are sampled are called ‘A-scans’, and the other scans 
are call ‘B-scans’. During the 102.4° arc of an A-scan, 64 samples of each of the lower 
channels are taken, with the integration period being 7.95 ms. The spatial resolutions of 
the samples depends upon the frequency and are given in Table 1. The IF passband for the 
19 and 22 GHz channels is 240 MHz. The passband for the 37 GHz channels is 900 MHz, 
and for the 85 GHz channel it is 1400 MHz. 
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Table 1 . Temporal and Spatial Resolution of SSM/I Channels 


I 

I 

I 


Frequency 

(GHz) 

Polarization 

Integration 

Period 

3 dB Footprint Size 
Along-track Cross-track 

19.35 

vertical 

7.95 ms 

69 km 

43 km 

19.35 

horizontal 

7.95 ms 

69 km 

43 km 

22.235 

vertical 

7.95 ms 

50 km 

40 km 

37.0 

vertical 

7.95 ms 

37 km 

28 km 

37.0 

horizontal 

7.95 ms 

37 km 

29 km 

85.5 

vertical 

3.89 ms 

15 km 

13 km 

85.5 

horizontal 

3.89 ms 

15 km 

13 km 












SECTION 3. SSM/I DEF TAPES 


We obtain the SSM/I data on Date Exchange Format (DEF) tapes. These tapes are 
originally produced at the Fleet Numerical Oceanography Center (FNOC), and copies are 
mailed to the National Environmental Satellite, Data, and Information Service (NESDIS). 
NESDIS then mails us copies of these DEF tapes. There are between 3 and 4 DEF tapes 
generated each day for each satellite. The large number of tapes is due to an inefficient 
format. We ingest the daily supply of DEF tapes, organize the data into a chronological, 
orbit-by-orbit data set, and output the data in a User friendly format onto 6250-bpi tapes. 
There are 8 T A tapes per month per satellite. We do not modify the DEF values for an- 
tenna temperatures or calibration data. However, we found problems with the cell latitudes 
and longitudes that are on the DEF tapes. Accordingly, starting with the January 1, 1989 
T a tape, we compute the cell latitudes and longitudes from a smoothed orbit ephemeris 
rather than using the DEF latitudes and longitudes. Prior to 1989, the T A tapes contain the 
DEF locations. The subroutine DECODE does the necessary location correction for these 
earlier tapes. Section 6 discusses the problems with the DEF locations and our methods for 
correcting these problems. 
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SECTION 4. COMPUTATION OF ANTENNA TEMPERATURES 


In computing the SSM/I antenna temperatures, the basic assumption is that the 
radiometer output voltage is linearly related to the input power at the mixer/preamplifier. 
Nonlinear effects such as imperfections in the square- law detector and IF amplifier com- 
pression are assume negligible. Expressing the input power in terms of radiation tempera- 
ture gives 


Vj - ck/lgO'^ + T n ) (1) 

where V i is the radiometer output voltage, (Kelvin) is the radiation temperature enter- 
ing the feedhorn, and T N (Kelvin) is the radiometer noise temperature. The constants are c 
= detector constant, k = Boltzmann constant, = receiver predetection bandwidth, and g = 
receiver gain. The subscript i denotes either a cold-space observation (i = C), a hot-load 
observation (i = H), or an Earth observation (i = E). The gain and noise temperature are 
assumed stable over the period of one scan, which is 1.9 s. That is to say, g and T N are as- 
sumed to be the same for the cold-space, hot-load, and Earth observations. The output 
voltage is converted to counts by applying a scaling factor and a constant offset: 


Cj = wV i -P (2) 

Combining (1) and (2) for i = C, H, and E gives 

Tae = A C e + B (3) 

A = (Tah " T AC )/(C h _ C c ) (4) 

B = (T AC C H " T ahCc)/(C h - C c ) (5) 


Henceforth, the subscript E on the Earth antenna temperature is dropped, and it is under- 
stood that T a refers to the Earth observation. For each SSM/I scan, five cold counts and 
five hot counts are recorded for each channel. (Only 85 GHz observations are taken on 
the B-scans). The quantities C c and C H are found by averaging the five individual counts. 
The cold space antenna temperature T AC is set to 2.7 K for all channels. The hot reference 
antenna temperature is found from 

Tah = Tqh + 0.01(T o p - T oH ) (6) 

where T oH and T oP are the physical temperatures of the hot load and the radiator plate 
facing the hot load. The small correction for the radiator plate is to account for radiative 
coupling between the hot load and the top plate of the rotating drum assembly which faces 
the hot load. There are three thermistors on the hot load and one on the radiator plate. 
The temperature T oH is found by averaging the three thermistor readings. The T A tapes 
contain the values for t a> T Q p, T oH , C c , and C H . Thus the A and B coefficients can be 
computed, and equation (3) can be inverted to recover the original values for the Earth- 
viewing counts C E . 

Starting with the F08 orbit 17057 on October 9, 1990, the FNOC procedure for com- 
puting T a was slightly modified. Rather than averaging the cold and hot counts for a 
single scan, the quantities C c and C H are found by averaging the counts for the last 20 
scans (10 A-scans and 10 B-scans). The purpose of this modification is to reduce the noise 
in C c and C H . Note that the last 20 scans are used regardless of data gaps. Thus, when 
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data gaps occur, the time period for this back-averaging will exceed 38 s, and some 
degradation in the calibration may occur. Also note that the number of scans in the back- 
average is less than 20 for the first 19 scans in a FNOC DEF file. We need to determine if 
this change in the calibration procedure produces a discontinuity in the F08 T A interannual 
time series. Any discontinuity will certainly be small (« 0.1 K), but climate change studies 
could still be affected. All F10 data have been processed with the new calibration proce- 
dure, and thus there is no discontinuity problem. 
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SECTION 5. COMPUTATION OF BRIGHTNESS TEMPERATURES 


Antenna temperature is a measure of the radiant power entering the feedhorn. It is 
computed by integrating over the gain pattern of the parabolic reflector and feedhorn as- 
sembly. 


T Ap ~ J (Gyp Tgy + G hp T Bh ) (7) 

where dn is the differential solid angle and the integral is over the entire 4tt steradians of 
a sphere. The subscript p equals h or v and denotes the feedhorn port. The terms and 
G h p are the antenna gains in the direction dO for vertically and horizontally polarized 
radiation received by port p. The terms T^ and T Bh are the vertically and horizontally 
polarized brightness temperatures coming from the dfl direction. The polarization is 
referenced to the Earth’s surface. Note that the gains G^ and G hp account for both the 
directional dependence of the antenna pattern and the rotation of the Earth- referenced 
polarization vector relative to the feedhorn polarization vector. 

Before launch, two properties of the antenna gain patterns for each SSM/I channel are 
measured: 1) the relative antenna gain (both co-pol and cross-pol) in the near-boresight 

direction and 2) the feedhorn spillover. The near-boresight gains are measured out to 6° 
for the 19 and 22 GHz channels, out to 4° for the 37 GHz channels, and out to 3° for the 
85 GHz channels. The feedhorn spillover is a measure of the power that enters the 
feedhorn directly from space, as opposed to the primary component of power that enters 
the feedhorn from the parabolic reflector. The feedhorn spillover is computed by measur- 
ing the percentage of the feedhorn gain pattern that is not subtended by the parabolic 
reflector. It is typically a few percent. Since the brightness temperature of space is much 
less than the Earth’s brightness temperature (2.7 K compared to 150-280 K), the spillover 
has an appreciable effect. 

In order to specify G^ and G hp for all directions, we model the antenna gain using the 
vector Kirchhoff approximation for a circular aperture illuminated by a plane wave, with 
an amplitude that tapers off according to the SSM/I feedhorn pattern. The effective radius 
of the aperture is chosen so as to match the SSM/I near-boresight measurements. Thus, the 
model agrees with the near-boresight and spillover measurements and provide the means to 
specify the antenna sidelobes that were not measured. 

The integral (7) is computed using an ocean brightness temperature model [Wentz, 
1991] to compute T^ and T Bh over the entire field of view as seen by the SSM/I at an al- 
titude of 860 km. The ocean environment is assumed uniform over the field of view. For 
the part of integral that corresponds to cold space, T^ and T Bh are set to 2.7 K. In this 
way, the antenna temperature is computed for various environmental conditions ranging 
from clear skies to light rain. An excellent approximation for all environmental conditions 
is found to be 


T Ap “ Qvp T Bvb + Qhp T Bhb + 2 - 7 Qop ( 8 ) 

where T^ and T Bhb are the v-pol and h-pol brightness temperatures in the antenna 
boresight direction (subscript b denotes boresight) and the Q terms are functions of the 
spillover factor 5 and the cross-polarization factor x p , P = v or h: 
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Q w = (1 - 6)/(l + X v ) (9) 

Qhv= X v (l -«)/(!+ Xv) (10) 

Q hh = (1 - 6)/(l + x h ) (11) 

Qvh - XhO - «)/(! + Xh) (12) 

Qov = Qoh = ^ (13) 

Note that the Q’s are normalized such that 

Qvp + Qhp + Qop = 1 (14) 

The spillover factor S is given by the integral of the antenna pattern over cold space 

5 = J dfl (G^ + G hp ) (15) 

space 


The leakage factor x p is a measure of the amount of radiation entering port p that has a 
polarization orthogonal to the port-p polarization mode. It is determined by fitting the 
above approximation (8) for T A to the exact integral computation of T A . Table 2 gives the 
values of 6 and x p for the seven SSM/I channels. For a given frequency, S is independent 
of the polarization. Over the considered range of environmental conditions, the error be- 
tween the T a approximation and integral is about 0.1 K. Note that the 0.1 K figure only 
indicates the accuracy of the linear approximation to the integral, assuming the modeled 
antenna patterns and a uniform ocean field of view. This does not mean that we know the 
antenna pattern accurately enough to compute T A to an absolute accuracy of 0.1 K. 

The S term is computed directly from the antenna patterns and does not depend on the 
assumed environmental scene. The x p term does have a small dependence on the choice of 
scene. For the range of ocean scenes considered (clear skies to light rain), the variation of 
Xp translates to about a 0.1 K variation in T B . Although the value of x p for land and ice 
observations differs slightly from the average ocean value given in Table P 2, this has negli- 
gible effect on the computation of brightness temperature because the T B polarization dif- 
ference for ice and land is small compared to water. As the scene becomes unpolarized, the 
T b computation becomes independent of x p > as is shown by equation (27) below. 


Table 2. 

Spillover and Polarization Leakage Factors 


Frequency (GHz) 

Polarization 

Spillover, 6 

Leakage, x p 

19.35 

vertical 

0.03199 

0.00379 

19.35 

horizontal 

0.03199 

0.00525 

22.235 

vertical 

0.02685 

0.00983 

37.0 

vertical 

0.01434 

0.02136 

37.0 

horizontal 

0.01434 

0.02664 

85.5 

vertical 

0.01186 

0.01387 

85.5 

horizontal 

0.01186 

0.01967 
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Before converting the antenna temperatures to brightness temperatures, we make cor- 
rections for two types of systematic errors. The first error seems to be due to the feedhorn 
partially seeing the cold-space reflector at the edge of the Earth- viewing portion of the 
scan. As a result, there is a systematic rolloff of the T A values starting near scan position 
100 and reaching a maximum value of about 1 K at the final scan position 128. This error 
was found by stratifying the SSM/I data according to scan position and then averaging over 
one year. The results for 1987, 1988, and 1989 are very similar. 

The second systematic error is the mean difference between the SSM/I antenna tem- 
perature observations and the Wentz [1991] radiative transfer model for the ocean and inter- 
vening atmosphere. The computation of this error is described in Wentz [1991]. The off- 
sets in terms of the observations minus the model are 2.0, 3.5, 1.3, -1.6, and -0.2 K for the 
19V, 19H, 22V, 37V, and 37H channels, respectively. These offsets are probably due to a 
combination of several errors including the following: 

1. Errors in specifying antenna parameters S and x p 

2. Errors in specifying the emissivity of the hot load 

3. Errors in the model formulation for the dielectric constant of sea water 

4. Errors in the model formulation for oxygen and water vapor absorption. 

Assuming that the along-scan error is due to the feedhorn seeing the cold-space reflec- 
tor, then the effect can be modeled simply as an error in the spillover term 1 - S. Ac- 
cording to T a equation (8), an error in 1 - 5 produces a multiplicative error in T A . It is 
less clear how to model the offset error because it is probably due to a combination of 
many effects. Accordingly, we simply assume that it is an additive error, and the following 
correction is applied ,to the antenna temperatures. 

TAcor = *1 Ta & (16) 

where T A is the uncorrected value computed from equation (3), T Acor is the corrected 
value, r) is the along-scan correction coefficient, and p is the offset listed in the previous 
paragraph. Figure 2 shows i? plotted versus scan position for the 37V channel for the F08 
satellite. The 19V, 19H, 22V, and 37H channels show a similar shape for the r? curve. The 
along-scan correction coefficient is normalized such that the value of r? average over all 
scan positions is equal to unity. 

Subroutine DECODE performs the along-scan correction if the User has selected the 
option for computing brightness temperatures (i.e., ITB = 1). However, DECODE does not 
perform the offset correction. That is to say, p is set to zero in DECODE. Section 11 
describes a simple modification to DECODE that will implement the offset correction. Ap- 
plying the offset matches the observations to the Wentz [1991] ocean model. However, this 
may not be appropriate for land and ice observations. 

No error analysis has been done for the 85 GHz channels, and no 85 GHz correction is 
applied. Also, the error analysis for the F10 satellite has yet to be done, and the current 
version of DECODE sets r) = 1 and p = 0 for F10. When the F10 error analysis is com- 
pleted, we will send the Users the tj and p values. 
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ALONG-SCAN ANTENNA TEMPERATURE CORRECTION COEFFICIENT 



Fig. 2. The F08 along-scan T A correction coefficient for the 37 GHz v-pol channel. 
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After correcting the antenna temperatures, the brightness temperatures are found. At 
19, 37 and 85 GHz, dual-polarization antenna temperature measurements are taken. 

"^Av = Qw "^Bv + Qhv "^Bh + 2,7 Q ov (17) 

T Ah = Qhh T Bh + Qvh T Bv + 2 - 7 Qoh U 8 ) 

where we have dropped the subscript b and it is understood that the brightness tempera- 
tures are in the boresight direction. It is also understood that the antenna temperatures T Av 
and are the corrected values given by equation (16). This system of two linear equa- 
tion is easily inverted to express brightness temperature in terms of the antenna tempera- 
tures: 



"^Bv A w Tav + A hv "^Ah + 2,7 A ov 

(19) 


^Bh — A hh "^Ah * A vh "^"av * 2,7 A oh 

(20) 

where the factor of 2.7 
of the spillover factor 5 

is the temperature of cold space. The A coefficients are 
and the cross-polarization factor x^ 

functions 


A w = (1 + X v )/[(1 - XyXhXl “ «)] 

(21) 


A hv = -X v 0 + Xh)/[d " X v X h )(l " 5)] 

(22) 


A hh ■ 0 + X h )/[(1 - X V Xh)0 - 5)] 

(23) 


A vh “ "XhO + X v )/[(1 - X v X h Xl - *)] 

(24) 


A ov - A 0 h = -«/U - 8) 

(25) 


At 22 GHz, only v-pol antenna temperatures are available. In order to derive an ex- 
pression giving T B as a function of T A at 22 GHz, we use the ocean brightness temperature 
model and the Q coefficients to compute T A for a wide range of environmental conditions. 
A least-squares regression yields the following relationship: 

Tbvt “ 1.01993 T Av + 1.994 (26) 

This proved to be a very accurate fit, with the rms error between the regressed T^ value 
and the actual value being 0.1 K, or less, for the range of antenna temperatures from 176 
to 270 K. The reason for the good fit is that the small cross-polarization h-pol leakage is 
highly correlated with the v-pol radiation. Note that this expression was derived for ocean 
observations. Over land, the brightness temperature is depolarized. Assuming that Tgy = 
T Bh , one can directly invert the T A equation, and obtain an alternative expression for T^ 
at 22 GHz. 

Tfiv = T Av /(1 - 6) - 2.7 5/(1 - 5) 


= 1.02759 T Av - 0.074 (27) 

Fortunately, these two expressions yield similar values in the 240 K. to 300 K range, where 
they differ by 0.2 K, or less. They intersect at 270 K. . Thus, the ocean-derived expression 
appears to be quite adequate for over land observations. A small error will occur for cases 
in which the brightness temperatures have relatively low values (below 240 K) and the 
polarization ratio is markedly different from the ocean polarization ratio. 
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The above computation of brightness temperature assumes a horizontally- uniform 
ocean scene. Any horizontal variability will be spatially smoothed by the antenna pattern, 
and as a result the T B computed from the above expressions will be a smoothed representa- 
tion of the actual brightness temperature field. Image enhancement techniques can be ap- 
plied to T b images to partially restore the actual T B contrast. The selection of the ap- 
propriate enhancement depends upon the particular application, and we leave this problem 
to the User. 
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SECTION 6. COMPUTATION OF 
ORBIT PARAMETERS AND EARTH LOCATIONS 


The cell latitudes and longitudes on the DEF tapes produced by FNOC are sometimes 
in error due to the following problems: 

1. There are some algorithm errors in the FNOC data processing software. 

2. The satellite ephemeris is sometimes incorrect due to spacecraft tracking errors and orbit 
prediction errors. This problem is particularly severe during periods of increased solar ac- 
tivity. 

3. The boresight nadir angle and the alignment of the SSM/I instrument relative to the 
spacecraft are slightly misspecified. 

There have been a number of FNOC algorithm errors. Sometimes the wrong Julian 
day is assigned to the orbit (see Section 7). Occasionally, the time for the spacecraft loca- 
tion and the scan time are very different. The two times should be within 60 s of each 
other, but at the begin and end of an FNOC DEF file, the two times can be as much as 
40,000 s apart. The spacecraft longitude is sometimes out-of-bounds. There are also other 
problems that have recently been corrected. For example, the ephemeris was not updated 
frequently enough, and the location of cell 128 was in error. 

In response to these problems, we developed our own routine for computing the 
latitudes and longitudes rather than using the DEF values. The input to our geolocation 
routine is the satellite ephemeris for a ±7 day period centered on the orbit being processed. 
The ephemeris is first subjected to quality control to removed any data that are clearly er- 
roneous. Next the ephemeris is smoothed to remove the noise. The smoothed ephemeris is 
then used to computed the SSM/I cell latitudes and longitudes. 

This procedure of computing latitude and longitude from a smoothed ephemeris was 
initiated with the January 1, 1989 F08 T A tape. Prior to 1989, the latitudes and longitudes 
on the T A tapes are the same as the values on the DEF tapes. Fortunately, an analysis of 
the ephemeris prior to 1989 does not show any large errors except for a brief period during 
the second half of February 1988. For these pre-1989 tapes, subroutine DECODE does the 
sensor alignment and the cell 128 correction while reading the tape. 

Our geolocation routine is free of the algorithm errors listed above except for the 
Julian day error. As each orbit is processed, we perform a number of cross-checks to 
verify the Julian day. Usually these cross-checks detect the Julian day error, which can 
then be corrected. However, a few of the Julian day errors will go undetected, and these 
‘erroneous’ orbits will be written to the T A tape. These erroneous orbits are later detected 
by the quality control routine discussed in Section 7. 

The third geolocation problem listed above is the misalignment of the SSM/I sensor 
relative to the spacecraft and the misspecification of the boresight nadir angle. For F10, 
this problem is corrected by our geolocation routine, and the locations on the T A tapes are 
correct. For F08, part of the correction is done by our geolocation routine and part is done 
by DECODE, as is explained in the next three paragraphs. 
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The F08 SSM/I requires a yaw adjustment and an along-track adjustment. The yaw 
correction is done by DECODE. The along-track correction is also done by DECODE for 
the T a tapes prior to 1989. For 1989 and thereafter, the along-track correction is done by 
our geolocation routine. Hence for the 1989 F08 T A tapes and thereafter, DECODE only 
does the yaw correction. DECODE automatically does the necessary corrections based on 
the scan time and the setting of the User option IADJ (See Section 11). 

The F08 yaw adjustment is a 0.5° clockwise rotation (looking down on the satellite) 
about nadir. For F08 T A tapes prior to March 1990, the along-track adjustment is a 15-km 
shift in the direction opposite to the spacecraft velocity vector. Thereafter, the adjustment 
is a 10-km shift in the same direction. The reason for the change is discussed in the next 
paragraph. These adjustments were derived by comparing SSM/I images with coastlines 
and are in good agreement with the corrections found by Hollinger [private communication, 
1991]. We used 85H images of the South Coast of Australia to determine the along-track 
error. The residual error after performing these corrections appears to be less than 10 km. 

In June 1989, a change was made to the ephemeris computation at FNOC. As a result 
the original 15-km along-track error was reduced to 10 km. For the T A tapes prior to 
1989, the 15-km adjustment is applied by subroutine DECODE if the User selects the op- 
tion IADJ = 1. For the T A tapes from January 1989 through February 1990, our geoloca- 
tion routine applies a 15 km translation when computing latitudes and longitudes. Starting 
with the March 1, 1990 T A tape, the geolocation routine applies a 10-km translation rather 
than a 15-km translation in order to be consistent with the FNOC June 1989 ephemeris up- 
date. Unfortunately there is a 9 month lag between the FNOC ephemeris update and our 
T a tape update. (We did not know that the new ephemeris changed the along-track error.) 
Thus from June 1989 through February 1990, there is a residual 5-km along-track error in 
the locations on the F08 T A tapes. 

For F10, the necessary adjustments are done by the geolocation routine, and the loca- 
tions on the T a tapes are correct. Two adjustments are done. First, the nadir boresight 
angle (see equation (45) below) is increased from its prelaunch value of 44.63° to 45.37°. 
Second, the azimuth start angle (see equation (46) below) is change from its prelaunch 
value of 50.9° to 51.1°. The change in the azimuth start angle is equivalent to a 0.2° 
counter-clockwise yaw adjustment. These F10 adjustments were derived by Hollinger 
[private communication, 1991]. 

In the remainder of this section, we give the formulation for computing the orbit 
parameters and Earth locations. We use a simple model for the SSM/I orbit. The 
spacecraft position at time t is specify by the vector R that extends from the Earth’s center 
to the spacecraft. This vector is given by its magnitude R and its three directional com- 
ponents R*, Ry, and Rj in a celestial coordinate system that does not rotate with the Earth. 
The x-axis is defined by the vector pointing from the Earth’s center to the celestial position 
corresponding to the most recent ascending node of the satellite (i.e., the point where the 
satellite last crossed the equator in the south-to-north direction). The z-axis is defined by 
the vector pointing from the Earth’s center to the North Pole. The y-axis is orthogonal to 
the x-axis and z-axis. The magnitude R (km) at time t is approximated by 

R = a [1 - € cos(T - T 0 )] + 1.82 (1 -2 sin 2 a sin 2 ^) (28) 

4> = T + 2 e [sin(T - T 0 ) + sin TJ (29) 

T = (2jr/P)(t - t 0 ) (30) 
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where a is the major semi-axis (km), e is the eccentricity, T 0 is the angle of perigee, P is 
the orbit period, a is the inclination angle, and t 0 is the time of the last ascending node. 
The second term in (28) with the factor 1.82 km is an empirical adjustment to account for 
the higher order orbit moments. Except for the empirical adjustment, the above expres- 
sions are simply Kepler’s equations expanded to first order in e. The eccentricity of the 
SSM/I’s orbit is less than 0.01, and the accuracy of this simple model is better than 1 km. 
In terms of the orbit angle <f>, the three directional components are 

R x = cos <f> (31) 

R y = cos a sin <f> (32) 

R z = sin a sin <p (33) 


The geocentric latitude 7 of R is 

7 = arctan[R z /(R x 2 + R y 2 ) 1/2 ] (34) 

The azimuth angle rp of R is defined as the angle between the x-axis and the projection of 
R onto the x-y plane. 

rp = arctan2(R y , R x ) (35) 

The spacecraft geodetic latitude and east longitude (in radians) are found from the follow- 
ing: 


s/c lat = arctan[(p e /p p ) 2 tan 7] (36) 

s/c Ion = ip + 2xft L - mod(t,86400)]/86400 (37) 

where p e and p p are the Earth’s equatorial and polar radii, which equal 6378.5 and 6356.9 
km, respectively. The parameter t L is the local time of the ascending node associated with 


the x-axis. The second term in (37) specifies the iongitude of the : 
Earth rotation. The satellite altitude is simply given by 

s-axis and accounts for 

s/c alt = R - p 

(38) 

P = PePp/[(Pp C0S i) 2 + (P e sin '7 ) 2 ] Vz 

(39) 

where p is the Earth’s radius along R. 



The DEF tapes give the spacecraft geodetic latitude, longitude, and altitude at ap- 
proximately 60 s intervals. For each orbit, the six orbit parameters (t 0 , a, t L , at, e, and T 0 ) 
are found by fitting (in a least-squares sense) the above expressions to the DEF latitudes, 
longitudes, and altitudes. Typical between 1% and 5% of the orbits have a poor quality of 
fit and are excluded from further processing. The parameters from the remaining orbits 
are processes by a regression routine that finds either a first-order or a second-order fit 
versus time to the parameter. First-order fits are used for parameters a, t L , and at, and 
second-order fits are used for the other three parameters. The fits are done over a ±100 
orbit window centered on the specified orbit. The resulting regressions are used to com- 
pute smooth values of the orbit parameters for each orbit, which are stored in an orbit 
parameter table. 
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Figures 3 and 4 show the smoothed orbit parameters for the F08 and F10 satellites, 
respectively. Instead of showing the ascending node time t 0 , we show the period P, which 
is the first derivative of t 0 . The figures show the variation of the parameters over 2000 or- 
bits that occurred during the first half of 1991. The dots are the values for the individual 
orbits that have a good quality of fit, and the solid curves are the regressions. 

The period P and the major semi-axis a show a slow, steady decrease with time. For 
example, at the beginning of the F08 mission in July 1987 the values for P and a were 
6118.5 s and 7229.3 km. Four years later, the value are 6111.0 s and 7223.3 km, respec- 
tively. The ascending node local time t L steadily increases with time. For F08 the change 
is small, being about 70 s/year. The F10 ascending node local time is increasing at a much 
faster rate of 2800 s/year. The inclination angle a is slowly decreasing with time at a rate 
of about -0.01 "/year for F08 and -0.05 8 /year for F10. The eccentricity e and angle of 
perigee T 0 are periodic. The period is about 1800 orbits and the amplitude remains con- 
stant over time. F10 has a much larger eccentricity than does F08 due to a malfunction 
that occurred during the F10 launch. (The period and the major semi-axis are redundant 
parameters in the sense that they are related by Kepler’s third law of motion. The expres- 
sion 21.6106 P 2 ^ 3 yields values for a that are within 0.2 km of the values in the orbit 
parameter table.) 

To compute the latitudes and longitudes along the SSM/I scan, it is first necessary to 
specify the spacecraft coordinate system, which is defined by three unit vectors: Uj, U 2 , 
and U 3 . The first vector is the nadir vector. It points from the satellite nadir point on the 
Earth up toward the satellite. The satellite nadir point is defined such that Uj is orthogonal 
to the Earth’s surface. Because the Earth is an oblate spheroid, Uj does not coincide with 
R except at the equator. In terms of the x,y,z coordinate system Uj is given by 

Uj = [cos r!> cos(t + vp/ R), sin tp cos(7 + vp/R), sin(7 + vp/R)] (40) 

The angle v (radians) is the difference between the geodetic minus geocentric latitude at the 
satellite nadir point, and p is the Earth’s radius at the satellite nadir point. These two 
quantities depend on the nadir point latitude. For the SSM/I altitude, the difference be- 
tween the nadir point geocentric latitude and 7 is less than 0.03°, and 7 can be used in 
place of the nadir latitude to very accurately approximate v and p. 

y = 1/2 [1 - (pp/p e ) 2 ] sin 27 (41) 

The second spacecraft coordinate vector U 2 is the spacecraft bow vector that points 
approximately in the direction of the spacecraft velocity. U 2 is orthogonal to Uj and lies in 
the orbit plane. It is given by 


U 2 = Uj x M/lUj x M| (42) 

where M is the vector that is orthogonal to the orbit plane. 

M = [0, sin a, -cos a] (43) 

The third spacecraft coordinate vector is 

U 3 = Uj x U 2 (44) 
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The antenna boresight unit vector B is given by 

B = -cos sin e n (U 2 cos w + U 3 sin u>) (45) 

where 0 n is the boresight nadir angle and w is the boresight azimuth angle. 

w = w 0 - 0.01396(i-l) (46) 

where i denotes scan position (i = 1 to 128) and a> 0 is the azimuth start angle for the first 
scan position. The factor 0.01396 is in terms of radians and corresponds to 0.8°. The 
boresight nadir angles for F08 and F10 are 44.75° and 45.37°, respectively. The F08 and 
F10 azimuth angles w 0 are 230.9° and 51.1°, respectively. 

Given R and B, the point where the extended boresight intersects the Earth surface is 
uniquely determined. Let C denote the vector from the center of the Earth to the 


boresight-Earth intersection point. This vector is given by 

C = R + rB (47) 

where r is the range distance from the satellite to the cell and is given by 

r = R[b - (b 2 - c) 1 *] (48) 

b - t- R A " R y B y - (P e /P P ) 2 R Z B Z ]/{1 + [(p e /p p ) 2 - 1 ]B Z 2 } (49) 

c = (1 - (p e /R) 2 + [(p e /p p ) 2 - 1 ]R z 2 )/{ 1 + [(Pg/Pp ) 2 - 1 ]B Z 2 } (50) 

where B x , B y , and B z are the x, y, and z components of B. The geodetic latitudes and east 
longitudes oi the cells are given by 

cell lat = arctan[(p e /p p ) 2 C z /(C x 2 + C y 2 ) 1/z ] (51) 

cell Ion = arctan2(C y , C x ) + 27r[t L - mod(t,86400)]/86400 (52) 

where C x , C y , and C z are the x, y, and z components of C. 


Cell latitudes and longitudes are found for every eighth scan position, i.e., i = 1,9, ..., 
121, and also for positions 123, 127, and 128. These 19 pairs of latitudes and longitudes 
are stored on the T A tape. When reading the tape, subroutine DECODE interpolates to find 
the locations of the other cells. The reason for not storing all locations on the T A tapes is 
to reduce the storage size. The interpolation used by DECODE is fast and is accurate to 
better than 2 km. 

The F08 yaw correction done by DECODE is essentially the same as subtracting 0.5° 
from w 0 . The F08 10-km along-track correction is done by adding 1.5 s to the ascending 
node times t 0 . It should be noted that a 0.2° increase in 6 n or a 0.2° decrease in the 
satellite pitch would also shift the location in the aft direction by about 10 km. It is dif- 
ficult to determine which is the best type of correction: timing, boresight angle, or pitch. 
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SECTION 7. ERRONEOUS DATA PERIODS 


Occasionally the F08 SSM/I data are mislocated by thousands of kilometers. The 
problem was first revealed by global images of the SSM/I T A ’s that showed the continent of 
South America up where Greenland should have been. What appears to be happening is 
that a wrong time tag is applied to the data coming down from SSM/I. This wrong time is 
then used to locate the data. In some cases it is clear that the time tag is off by exactly one 
day. Another problem with the F08 SSM/I data is that the T A values are sometimes com- 
pletely unphysical, ranging from 0 K to 650 K. The cause for this is unclear, but we 
suspect telemetry errors. 

We have not yet done an analysis of the F10 data and do not know if these two 
problems occur for F10. 

Both the mislocation problem and the out-of-bounds (OOB) problem occur throughout 
a DEF data file. A DEF data file usually corresponds to approximately one orbit. When 
doing statistical analyses of the SSM/I data, it is very important to exclude these time 
periods of erroneous data. Otherwise, the statistics can be significantly distorted. 

In order to detect periods of erroneous data, we compare the SSM/I T A observations 
with a global 1° latitude by 1° longitude map of averaged T A ’s observed during the previ- 
ous semimonthly period. A large and persistent discrepancy between the individual T A ’s 
and the map reveals the erroneous DEF data files. 

Figure 5 shows the differences between the individual F08 T A ’s and the average T A 
map for the first half of September 1987. Each horizontal line corresponds to a Greenwich 
day, with the Julian day of year shown at the left. The horizontal scale is Greenwich hours 
in the day. The vertical scale in the figure is degrees Kelvin, with the distance between 
horizontal lines being 150 K. 

Each point in the figure is the absolute value of a 90-second average of the difference 
between the 19-GHz n-pol T A observations and the semimonthly map. The 19-GHz h-poi 
channel was used because it provides the greatest sensitivity to varying surface types (land, 
water, ice) and to the zonal water vapor signature over the world’s oceans. 

Under normal conditions, the points lie near the horizontal lines, with an occasional 
small outlier. Storms, coast lines, and varying sea ice can cause larger differences, but only 
for a short period of time. The periods of mislocated orbits and OOB T A ’s are shown as 
clusters of widely varying points. Figure 5 shows five periods of erroneous data on Julian 
days 246, 250, 252, 253, and 255. 

Note that there are small, periodic features shown in the figure. These correspond to 
the SSM/I passing over coastlines and ice edges. The SSM/I revisits a particular feature 
every 14 orbits, which is 0.9912 days. Hence, as the Julian day increases the feature shifts 
to the left by 0.2 hours/day. Due to their periodicity, small size, and brief duration, these 
features can be distinguished from the mislocated orbits and OOB T A ’s. 

Starting with the second half of May 1988, the procedure for generating the figures 
was modified. For the OOB problem, we sometimes found that the 19-GHz h-pol channel 
showed no anomalous behavior, while one or more of the other channels were erroneous. 
That is to say, the OOB problem does not necessarily affects all channels. Also, an rms T A 
difference, rather than an absolute difference, proved to be a better discriminator of er- 
roneous data. Accordingly, we updated the error detection algorithm. 
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For each scan, the new algorithm first computes the mean squared difference (msd) 
between the T A ’s and the semimonthly map. Five msd’s are computed, corresponding to 
the SSM/I five lower frequency channels (i.e., 19V, 19H, 22V, 37V, and 37H). The maxi- 
mum msd is then selected, and these maxima are averaged over 90 seconds (24 A-scans). 
Let Y denote the resulting 90-second average. Normally, the value of Y is about 20 K or 
less. Accordingly, for displaying anomalies, 20 K is subtracted from Y. The result, 
Y' = Y - 20, is then plotted as a single point in the figures. Note that if Y' is less than 0 
K, it is set to 0 K; and if Y' is greater than 130 K, it is set to 130 K. 

A complete set of figures showing erroneous data periods for F08 can be obtained 
from Remote Sensing Systems. In addition to revealing periods of erroneous data, the 
figures provide a complete inventory of the SSM/I data set. They can be used to find 
periods of missing data. These data gaps are shown by the absence of plotted points. 
There is one figure for each semimonthly period from July 1987 to the present. For ex- 
ample, the AAUG87 figure denotes the August 1 through 15, 1987 period, and BAUG87 
figure denotes the August 16 through 31 period. We plan to produce similar figures for the 
F10 satellite. 

We have visually inspected all of the F08 figures and have identified the begin and end 
times for the periods of erroneous data. A small ASCII data file called BADLOC contain- 
ing these times is included on the enclosed 3.5" inch diskette. A listing of this data file is 
given in Appendix A. BADLOC is continually being updated. The enclosed version covers 
the first four years (July 1987 through June 1991) of F08 operation. Each record in the 
file contains 6 integers that specify the begin time and end time for the period of erroneous 
data. For example the first record is 

1987 198 4.0 1987 198 5.0 

The period of erroneous data begins in 1987 on Julian day 198 at Greenwich hour 4.0. The 
erroneous data ends one hour later at hour 5.0. These time windows are somewhat larger 
than the clusters of data outliers shown in the figures in order to ensure that all the er- 
roneous data are excluded. 

The problem of erroneous data was particularly bad during the second half of 1988. 
Then, in January 1989 and thereafter the occurrence of erroneous data greatly reduced. 
After 1988, there are, on the average, about one or two erroneous data periods each month. 
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SECTION 8. TIME SERIES OF RADIOMETER PARAMETERS 


The operation of the SSM/I is monitored by producing time series of various 
radiometer parameters. For example. Figure 6 shows the thermistor readings for three 
hardware components on the F08 SSM/I: the hot load, the top of the SSM/I drum, and the 
internal temperature of the electronics. Each of these component temperatures is shown as 
a different color in Figure 6 and is plotted versus time from July 1987 through June 1991. 
As can be seen, the temperatures go through a significant heating cycle every winter due to 
increase solar illumination. There is some concern that this large heating cycle may ad- 
versely affect the performance of the radiometers. The operation of the SSM/I radiometers 
is monitored by producing time series of the radiometer temperature resolution, noise tem- 
perature, and gain. The remainder of this section gives the formulation for computing 
these parameters and presents figures of the parameters plotted versus time for the first 
four years (July 1987 through June 1991) of the F08 operation. This type of analysis has 
not yet been done for the F10 satellite. 

The rms temperature resolution of the SSM/I channels is computed as follows. The 
variation of the radiometer calibration counts C (C denotes either cold counts C c or hot 
counts C H ) for a single A-scan is first found: 

£T C = [EC 2 - (EC) 2 /5]/4 (53) 

where E represents the sum over the five counts taken during an A-scan. The count 
variance is converted to a temperature variance by the following expression: 

cr T « A 2 a c (54) 

where the A coefficient is given by (4). The quantity <r T is then averaged over a period of 
3 or 4 days, and the rms temperature resolution is then given by 

AT = <c 7 ^ (55) 

where < > denotes the 3 or 4 day average. The temperature resolution found from the hot 
counts is slightly larger than that for the cold counts because the variation of the hot load 
emission is greater than the variation of the cold space emission. 

Figures 7 and 8 show the temperature resolution for the seven F08 SSM/I channels as 
found from the cold counts and from the hot counts, respectively. Each color represents a 
different channel. The most striking feature in Figures 7 and 8 is the degradation of the 
85V and 85H channels. Figure 9 is the same as Figure 7 except that the vertical scale has 
been compressed. Figure 9 can be used to determine those time periods for which the 85V 
and 85H channels are usable. The hash marks in these figures denote the beginning of the 
month. The horizontal axis goes from July 1, 1987 to July 1, 1991. 

During the first heating cycle in December 1987, the 85V temperature resolution 
begins to degrade. After the heating cycle is over, the 85V channel continues to degrade. 
After April 1988, the noise in the 85V channel exceeds 20 K. Between May 1988 and 
January 1989, there are time windows in which the noise subsides. The possibility exists 
that useful information can be obtained from the 85V channel during these time windows. 
After January 1989, the noise becomes very large, and the channel is not usable. 
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Fig. 7. F08 SSM/I radiometer temperature resolutions derived from the cold counts. 
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The 85H channel shows a slight degradation in the temperature resolution during the 
first heating cycle, but then recovers. The second heating cycle in December 1988 seems to 
cause a small but permanent degradation in the 85H channel. Then, the third heating cycle 
in December 1989 seriously damages the 85H channel, and the noise increases to between 5 
and 10 K. Curiously, during the fourth heating cycle in December 1990, the 85H channel 
temporarily recovers, but after January 1991 the noise becomes very large, exceeding 20 K. 

The other channels also experience a small increase in AT during the heating cycles. 
Fortunately, the lower frequency channels are not permanently affected by the heating, and 
AT returns to its normal level after the heating cycle is over. 

The radiometer noise temperature is found from equations (1), (2), and (4) to be 

T n = (Cj + P)A - T^ (56) 

where the A coefficient is given by (4). Either the cold sky observations or the hot load 
observations can be used in (56) (i.e., i = C or H), with identical results. The offset P is 
1782, 1796, 2835, 2614, 3178, 4167, and 4229 for the 19V, 19H, 22V, 37V, 37H, 85V, and 
85H channels, respectively. The relationship between the noise temperature and the 
radiometer temperature resolution is 

AT = (T m + T N )/09r)* (57) 

where f3 is the effective bandwidth of the radiometer and r is the integration time. We 
have verified that the temperature resolution computed from (57) is in good agreement with 
that computed from (55). Figure 10 shows T N for the seven channels. As expected from 
equation (57), this figure is very similar to Figures 7 and 8. 

The system gain G is defined by the reciprocal of (4). 

G = (C H - C c )/(T ah - T ac ) (58) 

This gain is the product of two components: 


G = G r G agc (59) 

where G r is the gain of the radiometer front-end and receiver and G agc is an attenuation 
applied by the automatic gain control (AGC) system. The variation in the receiver gain G r 
is primarily due to changes in the temperature of the receiver. The purpose of the AGC is 
to compensate for the G r variation such that the counts stay within the range of the digital 
converter. There are 16 AGC setting, and each setting adds 0.5 dB attenuation. When 
SSM/I was turned on (for the second time) in July 1987, the AGC system found initial set- 
tings for all channels. We denote these initial settings by G 0 . Thereafter, G agc was 
changed in steps of 0.5 dB in order to keep the counts within bounds. 

G = G r G 0 AG (60) 

where AG is the AGC setting relative to the initial AGC setting. The term AG is found 
from the AGC readings given on the T A tapes. In order to monitor the variation in the 
receiver gain relative to its initial gain, we computed the quantity G r G 0 by first computing 
G from (58) and then divided by AG. 
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Fig. 10. F08 SSM/I radiometer noise temperatures. 







Figure 1 1 shows G r G 0 for the seven SSM/I channels. The 85V and 85H gains decrease 
due to the heating cycles. The gain of the 37V and 37H channels also decreases with in- 
creasing temperature. In contrast, the gain of the 22V channel increases with increasing 
temperature. The gain of the 19V and 19H channels show little dependence on tempera- 
ture. 


The dependence of G r on temperature is one of the primary reasons for doing an ex- 
ternal calibration each scan. The temperature variation experience during a single orbit will 
produce changes in G r . Unless frequent calibrations are performed, the changing gain will 
introduce errors in the computation of the antenna temperatures. We have verified that the 
derivative of G r with respect to temperature displayed in Figure 1 1 is the same as that ob- 
tained from an analysis of the gain versus temperature variation for individual orbits. 
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SECTION 9. ORGANIZATION OF ANTENNA TEMPERATURE TAPES 

The SSM/I data reside on a series of 6250-bpi, 2400-ft magnetic tapes. A separate 
series of tapes is produced for each satellite. The first file on the tape is called the header 
file and contains general information about the tape. The remaining files are called data 
files and contain the SSM/I data. Each data file corresponds to a single SSM/I orbit. The 
beginning of an orbit is defined as the ascending equator crossing of the spacecraft (i.e., 
south- to -north crossing). The data files are chronologically ordered, and the last data file 
ends with a double end-of-file mark to specify the logical end of tape. 

For each month and each satellite, eight T A tapes are produced, except for a non- 
leap-year February for which there are only seven tapes. Table 3 shows the time periods 
corresponding to each of the tapes. A single tape contains either 3 or 4 days of data, ex- 
cept for February during a leap year, in which the eighth tape contains only 1 day of data. 
Each tape begins at the first ascending equator crossing for the Greenwich day specified in 
Table 3. 


Table 3. 

The Eight SSM/I T A 

Tapes per Month 

Tape Number 

Time Period (Day of Month) 


Except February 

February 

i 

1 through 4 

1 through 4 

2 

5 through 8 

5 through 8 

3 

9 through 12 

9 through 12 

4 

13 through 15 

13 through 16 

5 

16 through 19 

17 through 20 

6 

20 through 23 

21 through 24 

7 

24 through 27 

25 through 28 

8 

28 through end of 

' month 29 (for leap year) 


The orbit period is approximately 102 minutes, and thus there are 14.1 orbits per day. 
A tape containing three (four) days of data will have either 42 or 43 (56 or 57) data files, 
if there are no data gaps. A four-day tape with no data gaps contains about 160 

megabytes. 

The header file contains 23192 bytes divided into 13 logical records, except for the F08 
tapes prior to July 1, 1988, which contain only 21408 bytes divided into 12 logical records 
(see below). This record length of 1784 bytes is chosen to match the record length used for 
the data files. 

The T a data files consist of a series of chronologically-ordered logical records. Each 
record corresponds to a pair of A and B scans. The length of the logical record is 1784 
bytes. The logical records are grouped into tape data blocks, which are separated by a 0.3 
inch interblock gap. Each data block, except the last block in a file, contains 28544 bytes, 
which is 16 logical records. The last block in a file can contain any multiple of 1784 bytes, 
but not exceeding 28544 bytes. The logical record format for the T A data is described in 
Section 10. 
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The remainder of this section discusses the contents of the header file. The first 12 
records in the header file each contain 1784 ASCII characters. The following is an example 
of these 12 ASCII records 

SSM/I F08 TAPE 1991_JUN_P2_A, 55 DATA FILES -91129320242.0 

SPEC. BEG DAY 156 Wed Jun 5, 1991 00:15:10 

ACTUAL BEG DAY 156 Wed Jun 5, 1991 00:15:15 

SPEC. END DAY 160 Sun Jun 9, 1991 01:00:43 

ACTUAL END DAY 160 Sun Jun 9, 1991 01:00:42 

SPECIFIED ORBITS: 020427 • 020483 

ACTUAL ORBITS: 020427 - 020483 

BEGIN GAP = .08 MIN, END GAP = .02 MIN 

PERCENTAGE OF MISSING DATA - 16.389 % 

PRODUCED BY REMOTE SENSING SYSTEMS 
1101 COLLEGE AVE., SANTA ROSA, CA 95404 
CONTACT FRANK WENTZ 707-545-2904 (F.WENTZ/0MNET) 

The 12 records are padded on the right side with trailing ASCII blanks. This example is 
the header file for the second F08 T A tape in June 1991. The first line identifies the satel- 
lite and the time period according to year, month, and part of month (i.e., PI through P8). 
The first line also gives the number of data files on the tape that follow the header file. 
The last number on the first line is a checksum. It is the sum of all 2-byte words in all 
data files on the tape. We use this number to verify tape copies. 

There have been two changes to the first line in the header file. For the F08 tapes 
prior to 1991_APR_P7, the satellite number is not specified, and the format for the first 
line is 

COMPACT TA TAPE 1991_APR_P6_A, 57 DATA FILES -133747836598.0 

The first header line for all F10 tapes specifies the satellite number. Thus if the satellite 
number is not specified, then the tape is an F08 tape. For the F08 tapes prior to 
1989_SEP_P1, the checksum is not included, and the format is 

COMPACT TA TAPE 1989_AUG_P8_A, 56 DATA FILES 

The next four lines in the header file give the specified and actual begin and end times 
for the tape. The Julian day of year is given along with the date and time. The specified 
begin time is the time of the first ascending equator crossing specified for the T A tape. 
The specified end time is the time of the first ascending equator crossing specified for the 
next T a tape. The actual begin and end times are the first and last SSM/I scan times on 
the tape. Because of data gaps, the actual begin time may occur later than the specified 
begin time, and the actual end time may occur sooner than the specified end time. The 
next two lines in the header file give the specified and actual range of orbits. The dif- 
ferences (i.e., time gaps) between the specified and actual begin and end times are then 
given, and the following line shows the percentage of missing data for the specified time 
period. The remaining lines give the name and address of the organization responsible for 
producing the tape. 

The 13th logical record is an inventory of the data files. Note that this inventory is 
not on the F08 tapes prior to 1988_JUL_P1. For these earlier tapes, there are only 12 
records in the header file. The inventory record lists each orbit on the tape along with the 
times for the first and last scans in the orbit. 
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The following Fortran code can be used to print out the header file along with the tape 
inventory information. 


SUBROUTINE HEADER (LREC) 

LREC IS THE ARRAY THAT CONTAINS THE BYTES IN THE TAPE BLOCK 
THE USER MUST SUPPLY THE TAPE READ ROUTINE FOR FILLING LREC 
BEFORE CALLING THIS ROUTINE 

CHARACTERS LREC(1784,16> 

INTEGERS IBUF(1784) 

DATA N1,N2,N3/256,65536, 16777216/ 

PRINT OUT FIRST 12 LINES OF HEADER FILE 

DO 100 IREC=1,12 

URITE(6,6001) (LREC( I , IREC), 1=1 , 72) 

>001 FORMAT(1X,72A1) 

100 CONTINUE 

PRINT OUT INVENTORY FOR TAPE 

DO 200 1=1,1784 
IBUF(I)=ICHAR(LREC(I,13>) 

200 CONTINUE 

N=— 1 1 

DO 300 1=1,100 
N=N+12 

IREV IS THE ORBIT NUMBER 

ITIME1 IS THE TIME FOR THE FIRST SCAN IN ORBIT 

ITIME2 IS THE TIME AT THE LAST SCAN IN THE ORBIT 

IT1ME1 AND ITIME2 ARE IN TERMS OF SECONDS FROM BEGIN OF 1987 

IREV= N3*IBUF(N) ♦N2*IBUF(N+1)+N1*IBUF<N+ 2)+IBUF(N+ 3) 
ITIME1= N3*IBUF(N«-4)+N2*IBUF(N+5)+N1*IBUF(N+ 6)+IBUF(N+ 7) 
ITIME2= N3*IBUF<N+8)*N2*IBUF(N+9)+N1*IBUF(N+10)+IBUF(N+11) 

IREV=0 MEANS YOU HAVE COME TO THE END OF THE INVENTORY 

IF(IREV.EQ.O) GO TO 310 
URITE(6,6002) IREV, ITIME1, ITIME2 
6002 FORMAT(1X, 16,2111) 

300 CONTINUE 
310 CONTINUE 
RETURN 
END 
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SECTION 10. LOGICAL DATA RECORD FORMAT 


The information in a logical data record is stored as a string of 1784 bytes, with each 
byte representing an integer ranging from 0 to 255. The subroutine DECODE converts this 
byte string to SSM/I sensor parameters (i.e., spacecraft and cell locations, calibration data, 
and antenna temperatures). The sensor data are store in the common /OUTDAT/, which 
can be interfaced with the User’s application program. A description of the DECODE call- 
ing procedure and the contents of /OUTDAT/ is given in Section 11. 

We strongly recommend that DECODE be used to read the T A tapes, rather than the 
User developing his/her own routine. The User does not need to know the logical record 
format for the T A tapes because DECODE performs all the necessary conversions. The 
User only needs to call DECODE and then access the common /OUTDAT/. However, for 
some specialized applications, information on the record format may be useful. Accord- 
ingly, this section describes the data record format. 

Table 4 shows the byte format for a logical record. The record is divided into 29 
fields, with each field corresponding to a different set of data items. The table gives the 
number for the first byte in the field, the number of bytes in the field, the number of 
items in the field, the bias and scale factors to convert to the units specified, and a brief 
description of the data item. For example, Field 8 starts at byte 29, is 6 bytes long, and 
contains 3 hot load temperatures each stored as a 2-byte integer. 

For Fields 1 and 2 and Fields 4 through 26, the conversion from tape bytes to the ap- 
propriate sensor parameter is simple. For these fields, each item corresponds to either 2 or 
4 bytes, as indicated by Table 4. The bytes are first converted to a positive 16-bit or 32- 
bit integer I, where the first byte on the tape is the most significant. Letting S and B 
denote the scale and bias values given in Table 4, the sensor parameter P is given by 

P = S (I - B) (61) 

where P is in terms of the units specified in Table 4. 

The contents of Field 3 have undergone two changes. Originally, Field 3 was used to 
stored the time for the spacecraft location, i.e., the time associated with the spacecraft 
latitude, longitude, and altitude. Beginning in 1989, the spacecraft location is computed at 
the scan time given in Field 1, hence it was not necessary to have a separate time for the 
spacecraft location. Thus, starting with the January 1, 1989 T A tape, Field 3 was used to 
store an orbit number that was more precise than the value stored in Field 2. Then starting 
with the September 1, 1989 T A tape the precise orbit number was stored in Field 2, thereby 
freeing Field 3. Finally, starting with the August 1, 1989 T A tapes for F08 and F10, Field 
3 is used to stored the satellite identification number and the incidence angle for the scan. 

DECODE determines which T A tape is being processed by looking at the scan time and 
then correctly interprets of the contents of Field 3. Thus, the Field 3 format changes are 
transparent to the User. 

The contents of Fields 13 through 19 have also undergone a change. Prior to Septem- 
ber 1989, these fields contained various types of diagnostics data, which are of no use to 
the Users. Starting with the September 1, 1989 T A tape, these fields contain the SSM/I or- 
bit parameters listed in Table 4. For the T A tapes prior to September 1989, DECODE sets 
the seven orbit parameters to zero. 
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TABLE 4. DESCRIPTION OF LOGICAL DATA RECORDS ON SSM/I ANTENNA TEMPERATURE TAPES 
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Note 1: For T A tapes prior to January 1, 1989, Field 3 is the time for the spacecraft location. For T A tapes from January 1, 1989 through July 31, 
1991, Field 3 is the precise orbit number. For the August 1, 1991 T A tape and thereafter, Field 3 contains the satellite identification number and the 
incidence angle. 

Note 2: For T A tapes prior to September 1, 1989, Fields 13 through 19 contain diagnostic data. Thereafter, orbit parameters are stored in the fields. 


For Fields 27, 28, and 29, the data are stored in special packed formats, and the above 
conversion (61) is not applicable. Each pair of latitude and longitude differences in Field 
27 is stored as a single 2-byte integer, which must be decoded in a special way. The an- 
tenna temperatures in Fields 28 and 29 are stored as 12-bit integers, and the surface-type 
indices in Field 28 are stored as 3-bit integers. The procedure for unpacking these data is 
given in the listing of subroutine DECODE in Appendix B. 

The T a tapes do not contain locations for all of the A-scan and B-scan cells. Rather, 
19 pairs of latitudes and longitudes are stored for each scan. The subroutine DECODE in- 
terpolates these 19 locations to find the locations for the intermediate cells. The interpola- 
tion error is less than 2 km. 

More information on the logical record format is given in Section 11, which provides a 
cross-reference between the fields shown in Table 4 and the SSM/I sensor parameters that 
are stored in common /OUTDAT/. 
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SECTION 11. DESCRIPTION OF SUBROUTINE DECODE 


A Fortran subroutine named DECODE is provided to the User for decoding the infor- 
mation on the T a tapes. The subroutine is written in Fortran 77 and can be easily imple- 
mented on most computer systems. Before calling DECODE, the User must use his/her 
own routine to read from the tape. The bytes in a tape block must be read into the follow- 
ing array: 

CHARACTER* 1 LREC(1784,16) 

The first dimension corresponds to the number of bytes in a logical record, and the second 
dimension corresponds to the maximum number of logical records in a tape block. See Sec- 
tions 9 and 10 for further information on the tape format. Once LREC is filled, the User 
calls DECODE using the following statement: 

CALL DEC00E( I85GHZ, ITB, IADJ, IREC,LREC) 

The definitions of the User-defined arguments are given below. DECODE performs a 
number of functions. The bytes in LREC are converted to integers, which are then com- 
bined to compute values for the various SSM/I data items, including the spacecraft location, 
the calibration data, the cell locations, and the antenna temperatures. If ITB=1, then 
brightness temperatures are computed instead of antenna temperatures. These quantities are 
stored in a common area called /OUTDAT/, and program control is then returned to the 
User’s program. 

It is important for the User to understand the SSM/I sampling scheme. Each SSM/I 
scan is divided into 128 positions, with the distance between positions being 12.5 km. The 
along-track distance between scans is also 12.5 km. The 85 GHz channels are sampled at 
each scan position. Thus, the 85 GHz observations are available on a 12.5 km grid. A dif- 
ferent sampling scheme is used for the five lower channels (i.e., 19V, 19H, 22V, 37V, and 
37H). Observations of the lower channels are only taken every other scan. The scans 
during which the lower channels are sampled are called 5 A-scans’, and the scans during 
which only the 85 GHz channels are sampled are called ‘B-scans’ (see Figure 1). The dis- 
tance between A-scans is then 25 km. The integration time for the five lower channels is 
twice that of the 85 GHz channels, and hence each lower channel is only sampled 64 times 
over the A-scan. Thus, the five lower channels are available on a 25 km grid. The phase 
of the integration time for the five lower channels is such that the locations of the 64 
samples correspond to the locations of the odd 85 GHz scan positions. For example, the 
location of the 32 th sample of the lower channels is the same as the location of the 63 lh 
sample of the 85 GHz channels. Each logical record on the T A tape contains data for a 
single pair of A- and B-scans. 

The remainder of this section is divided into three part. The first part discusses the 
specification of the arguments for DECODE. The second part is an overview of the execu- 
tion of DECODE and the five subroutines called by DECODE. The third part is a detailed 
description of the common /OUTDAT/ 
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Subroutine DECODE has five arguments that must be specified by the User. The 
definition of these arguments follows: 

I85GHZ 

This argument must be set to either 0 or 1. The User should set I85GHZ=0 if information 
on the two 85 GHz channels is not required. When I85GHZ=0, DECODE does not process 
the 85 GHz data, thereby substantially reducing the processing time. In this case, DECODE 
does not compute values for the B-scan arrays BLAT, BLON, BTAHI, and IBTOIL. Also, 
DECODE does not find values for the A-scan array ATAHI and for the even-numbered 
elements in the A-scan arrays ALAT, ALON, and IATOIL. These arrays are defined 
below. The even-numbered elements on the A-scan correspond to those cells for which 
only 85 GHz observations are taken. When I85GHZ is set to 1, all channels are processed. 

ITB 

This argument must be set to either 0 or 1. When ITB=0, the arrays TALO, ATAHI, and 
BTAHI are filled with the antenna temperatures computed from equation (3). When ITB=1, 
these arrays are filled with brightness temperatures computed from equations (19), (20), and 
(26). Also when ITB=1, the along-scan correction r? specified by equation (16) is applied, 
but the offset correction p is set to zero. The offset correction removes small biases be- 
tween the T A ’s and the Wentz [1991] radiative transfer model for the ocean. Along-scan 
and offset T A corrections have not yet been found for the F10 SSM/I. If the User wants to 
implement the offset correction for F08, a simple change to the following two lines in 
subroutine FDTB08 is required. Comment out the first line and activate the second line. 

DATA ABIAS/5*0./ 

C DATA ABIAS/2. 0,3. 5, 1.32, -1.57,-0.22/ 

For F08, if the User wants the 85H brightness temperatures for time periods when the 85V 
antenna temperatures are in error (see Section 8 and Figure 9), then a special change in 
subroutine FDTB08 is required to prevent the 85V T A error for corrupting the 85H T B 
computation. Set the 85H cross-polarization leakage term .01967 to 0. in the following line. 

DATA CHI/. 00379, .00525, .00983,0.0, .02136, .02664, .01387, .01967/ 


IADJ 

This argument must be set to either 0 or 1. If IADJ=1, then DECODE does a 15-km 
along-track correction for the F08 tapes prior to 1989. No correction is done when 
IADJ=0. The setting of IADJ only affects the processing of F08 tapes prior to 1989. That 
is to say, for the F08 T A tapes after 1988 and for all F10 tapes, DECODE ignores the set- 
ting of IADJ because no along-track correction is needed for these tapes. Note that the 
along-track correction considerably slows down the processing speed, and hence IADJ 
should only be set to 1 if location accuracy is important. The correction translates the 
SSM/I imagery in the aft direction by 15 km. 

IREC 

This argument specifies which logical record is to be processed by DECODE. It must be 
set to a value between 1 and 16, and its value should not exceed the maximum number of 
logical records last read from the T A tape. Note that the last block in a tape file can con- 
tain less than 16 logical records. 

LREC 

This is the character array described above. It contains the bytes that the User’s program 
read from a data block on the T A tape. 
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We now give an overview of the execution of DECODE and the five subroutines called 
by DECODE. 

SUBROUTINE DECODEn85GHZ. ITB. IADJ. TREC. I.RF.Ct 

The arguments for this routine are described above. This routine begins by converting 
the specified logical record in the character array LREC to the integer array IBUF, which 
is used for all subsequent calculations. The routine then finds the spacecraft time, location, 
and orbit number. The satellite number ISAT (8=F08 and 10=F10) and incidence angle are 
next found. The calibration data and orbit parameters are then computed. DECODE next 
determines if an along-track adjustment is required. This adjustment is only required if 
the User has specified IADJ=1 and if the F08 scan time is before 1989 (i.e., IOLD=l). 
Note that the along-track correction requires the cell locations on both the A-scan and 
B-scan. Hence J85GHZ is set to 1 when the along-track correction is required. FDLTLN 
is then called to find the cell latitudes and longitudes. If the User has specified I85GHZ=1, 
then J85GHZ is set to 1 and FDLTLN finds locations for all cells on both scans. Other- 
wise, FDLTLN finds only the locations of the odd positions cells on the A-scan. The sub- 
routine FDTA is next called to compute the antenna temperatures. Finally, if the User has 
specified ITB=1, DECODE calls FDTB08 if the satellite is F08 or FDTBOO if the satellite is 
not F08 to convert the antenna temperatures to brightness temperatures. 

SUBROUTINE FDLTLNf J85GHZ. IOLD. YAW! 

This routine computes the latitudes and longitudes for the cells along the SSM/I scan. 
If the first argument J85GHZ=0, FDLTLN finds only the locations for the odd scan posi- 
tions on the A-scan. If J85GHZ=1, then locations are found for all scan positions on both 
the A-scan and the B-scan. If the second argument IOLD=l, then a correction is applied 
to the location of scan position 128. This correction is only done for the F08 tapes prior to 
1989. The third argument YAW specifies the yaw correction. This argument is scaled such 
that YAW=1 represents a 1.6° clockwise rotation (looking down on satellite) about nadir. 
The T a tapes do not contain locations for all cells. Rather, the latitudes and longitudes for 
19 scan base positions are stored on tape, and the locations for the other cells are found by 
means of an interpolation. The error in the interpolation is less than 2 km. FDLTLN first 
finds the latitudes and longitudes for the 19 A-scan base positions. The remaining A-scan 
locations are then found by interpolating. If J85GHZ=0, the yaw correction is applied to 
the A-scan, and control returns to DECODE. If J85GHZ=1, the 19 B-scan base positions 
are found. The T A tapes contain the 19 differences of the B-scan latitudes and longitudes 
minus the A-scan latitudes and longitudes. Storing the differences requires less bytes than 
storing the full B-scan latitudes and longitudes. The 19 differences are added to the 19 
A-scan base locations to obtain the 19 B-scan base locations. The interpolation is then 
done to obtain the locations at the other B-scan positions. The yaw correction is then ap- 
plied to both the A-scan and B-scan locations, and control returns to DECODE. 

SUBROUTINE ADJLOCITRK AD.I1 

This routine shifts the cell locations in a direction parallel to the spacecraft velocity vector. 
The argument TRKADJ is scaled such that a value of TRKADJ=1 (-1) corresponds to a 
12.5 km shift in the direction of (opposite to) the velocity vector. 
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SUBROUTINE FDT A(I85GHZI 

This routine computes the antenna temperatures. If the argument I85GHZ=0, then the 
85 GHz antenna temperatures are not computed. Also if I85GHZ=0, the routine finds only 
the surface-type indices for the odd A-scan cells. If I85GHZ=1, then all antenna tempera- 
tures and surface-type indices are found. The antenna temperatures are stored on the T A 
tapes as 12-bit words, and the surface-type indices are stored as 3-bit words. Most of the 
code in FDT A is for the unpacking of these bits. The precision of the T A ’s is 0.1 K except 
for anomalous values that exceed 380 K, for which the precision is 1.0 K. 

SUBROUTINE FDTB08n85GHZl 

This routine applies corrections to the F08 antenna temperatures and then converts the cor- 
rected antenna temperatures to brightness temperatures. If the argument I85GHz=0, then 
the 85 GHz channels are not processed. The first time the routine is called, it computes the 
T A -to-T B conversion coefficients and the along-scan correction arrays SBIAS1, SBIAS2, 
SBIAS3, SBIAS4, and SBIAS5. FDTB08 applies the along-scan correction rj and the offset 
correction p to the T A ’s, except that no correction is applied to the 85 GHz channels. The 
corrected T A ’s are then converted to T B ’s and control returns to DECODE. 

SUBROUTINE FDTB00CI85GHZ1 

This routine is the same as FDTB08 except that the T A correction arrays ABIAS, SBIAS1, 
SBIAS2, SBIAS3, SBIAS4, and SBIAS5 are set to zero. This routine is currently being used 
for the F10 and subsequent satellites until the T A corrections are found. 
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The last part of this section gives a detailed description of the contents of the common 
/OUTDAT/. The form of /OUTDAT/ is: 

c 

REAL*8 REV, XT I ME 

l NTEGER** I T I ME , l T I MSC , I VOLT , I AGC , l COLOA , I HOTA , I COLDB, IHOTB, I ASCTM 
INTEGER** 1ATOIL, IBTOll, ISAT 

REAL*4 XLATSC , XLONSC , ALTSC, THT , HLTEMP, RFTEMP , FRTEMP 
REAL** PERIOO,ASCLOC,ANGINC, AXIS, ECC.ANGPER, SPACER 
REAL** AL AT , ALON , BLAT , BLON , TALO, AT AH I , BTAH I 
C 

COMMON/OUTDAT/ REV, XTIME, ITIME, ITIMSC,XLATSC,XLONSC, ALTSC, THT, 

1 HLTEMP<3), IVOLT(2), RFTEMP, FRTEMP, IAGC<6), 

2 I ASCTM, PER IOO, ASCLOC,ANGINC,AXIS,ECC,ANGPER, ISAT, SPACER(6), 

2 IC0LDA(5,7), IH0TA(5,7), ICOLOB<5,2), IHOTB(5,2), 

3 ALAT < 1 28) , ALON ( 1 28 ) , BLAT < 1 28) , BLON< 1 28 ) , 

* TAL0(5,6*),ATAHI (2, 128),BTAHI(2, 128), IAT0IL(128), IBT0ILO28) 

The following is a description of each variable and array in /OUTDAT/. At the end of 
each description, we specify which field(s) in the Table 4 logical data record structure is 
used to compute the parameter. 

REV 

This 8 -byte real variable is the orbit position of the spacecraft. Prior to 1989, the time 
corresponding to REV is given by ITIMSC. For 1989 and thereafter, the time correspond- 
ing to REV is given by XTIME. REV equals an integer value when the spacecraft is 
crossing the equator in the south-to-north direction (i.e., the ascending node). For ex- 
ample, REV= 100.5 indicates that the spacecraft is halfway through orbit 100 and is crossing 
the equator in the descending, north-to-south direction. {Field 2 or Field 3} 

XTIME 

This 8-byte real variable is the Greenwich time at which the SSM/I B-scan begins. The 
time for the beginning of the A-scan is found by subtracting 1.9 s. This time is in terms 
of seconds from the beginning of 1987. Hence XTIME=0 denotes 00:00:00 Z, January 1, 
1987. For the F08 T A tapes prior to 1987_SEP_P6, the fractional part of XTIME was not 
computed, and hence XTIME is the nearest integer to the actual time. That is to say, 
XTIME=ITIME for the early F08 T A tapes. {Fields 1 and 5} 

ITIME 

This 4-byte integer is the nearest integer to XTIME, i.e., ITIME=NINT(XTIME). {Field 1} 
ITIMSC 

This 4-byte integer is the nearest integer to the Greenwich time corresponding to the 
spacecraft location as specified by the spacecraft latitude XLATSC, longitude XLONSC, 
and altitude ALTSC. For T A tapes prior to 1989, this time is usually within 60 seconds of 
the scan time XTIME. Starting with the January 1, 1989 T A tape and thereafter, the 
spacecraft location is calculated at the scan time XTIME, and hence ITIMSC equals ITIME. 
This time is in terms of seconds from the beginning of 1987. {Field 3 for 1987 and 1988; 
Field 1 thereafter) 
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XLATSC 

This 4-byte real variable is the spacecraft geodetic latitude in units of degrees. Negative 
latitudes correspond to south latitudes, and positive values correspond to north latitudes. 
This latitude corresponds to the location where a line from the Earth’s center to the satellite 
intersects the Earth’s surface. It is given by equation (36). Prior to 1989, the time cor- 
responding to XLATSC is given by ITIMSC. For 1989 and thereafter, the time correspond- 
ing to XLATSC is given by XTIME. {Field 4) 

XLONSC 

This 4-byte real variable is the spacecraft east longitude in units of degrees. This longitude 
corresponds to the location where a line from the Earth’s center to the satellite intersects 
the Earth’s surface. It is given by equation (37). Prior to 1989, the time corresponding to 
XLONSC is given by ITIMSC. For 1989 and thereafter, the time corresponding to 
XLONSC is given by XTIME. (Field 6} 

ALTSC 

This 4-byte real variable is the spacecraft altitude in units of kilometers. It is the distance 
from the satellite to the location where a line from the Earth’s center to the satellite inter- 
sects the Earth’s surface. It is given by equation (38). Prior to 1989, the time correspond- 
ing to ALTSC is given by ITIMSC. For 1989 and thereafter, the time corresponding to 
ALTSC is given by XTIME. {Field 7} 

THT 

This 4-byte real variable is the incidence angle made by the antenna boresight vector and 
the vector normal to the Earth’s surface at the point where the boresight vector intersects 
the Earth’s surface. It is in units of degrees. It is computed as a function of the spacecraft 
altitude ALTSC and the Earth’s radius of curvature at latitude XLATSC. The accuracy of 
the approximation being used is 0.03°. Prior to 1989, the time corresponding to THT is 
given by ITIMSC. For 1989 and thereafter, the time corresponding to THT is given by 
XTIME. {Fields 4 and 7 or Field 3} 

HLTEMPI31 

This 4-byte real array contains the temperatures recorded by the three thermistors attached 
to the SSM/I external hot load. These temperatures are in units of Kelvin degrees. The 
average of these three temperatures is denoted by T oH in equation (6). {Field 8} 

IVOLTm 

This 4-byte integer array contains two reference voltages for the SSM/I receiver. These 
voltages are in terms of counts and are not required by the calibration algorithm. {Field 9} 

RFTEMP 

This 4- byte real variable is the temperature recorded by the thermistor attached to the 
SSM/I R.F. mixer. It is in units of Kelvin degrees and is used to monitor the internal tem- 
perature of the receiver. It is not required by the calibration algorithm. {Field 10} 

FRTEMP 

This 4-byte real variable is the temperature recorded by the thermistor attached to the top 
plate of the rotating drum assembly which faces the hot load. It is in units of Kelvin de- 
grees. FRTEMP is denoted by T oP in equation (6). {Field 11} 
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IAGCI6) 

This 4-byte integer array contains the automatic gain control (AGC) readings in units of 
counts. The first three elements are the three reading for the A-scan, and the last three 
elements are the readings for the B-scan. For the A-scan, the AGC values for the in- 
dividual channels are given by the following expressions: 

IAGC85V = I NT ( I AGC( 1 )/256) 

IAGC85H = INT(( I AGCO ) -256*1 AGC85V)/16) 

IAGC37V = I AGCO )- 256*1 AGC85V- 16*IAGC85H 
IAGC37H = INT ( IAGC(2)/256) 

IAGC19V = INTO I AGC(2) -256*1 AGC37HJ/16) 

IAGC19H = I AGC ( 2 ) • 256* I AGC37H • 1 6* I AGC 1 9V 
IAGC22V = INTOAGC(3)/2S6) 

The B-scan values are also given by the above expressions, except that IAGC(4), IAGC(5), 
and IAGC(6) replace IAGC(l), IAGC(2), and IAGC(3), respectively. The AGC values are 
for diagnostic purposes only and are not required by the calibration algorithm. {Fields 12 
and 22} 

IASCTM 

This 4 -byte integer is the Greenwich time of the ascending node for the current orbit. The 
ascending node time is defined as the time at which the satellite crosses the equator in the 
south- to- north direction. It is denoted by t 0 in Section 6. This time is in terms of seconds 
from the beginning of 1987. For T A tapes prior to September 1, 1989, this orbit parameter 
is not available and is set to zero. {Field 13} 

PERIOD 

This 4-byte real variable is the period for the current orbit. The period is defined as the 
time between one ascending node and the next. It is denoted by P in Section 6. The 
period is in units of seconds. For T A tapes prior to September 1, 1989, this orbit parameter 
is not available and is set to zero. {Field 14} 

ASCLOC 

This 4-byte real variable is the local time of the ascending node for the current orbit. It is 
in terms of seconds in a day. It is denoted by t L in Section 6. The east longitude of the 
ascending node is given by 

XLONASC = 360. *( ASCLOC • MOO(IASCTM,86400))/86400. 

IF (XLONASC.LT .0. ) XLONASC=XLONASC+360. 

For T a tapes prior to September 1, 1989, this orbit parameter is not available and is set to 
zero. {Field 15} 

ANGINC 

This 4-byte real variable is the inclination angle for the current orbit. It is in units of de- 
grees and is denoted by a in Section 6. For T A tapes prior to September 1, 1989, this orbit 
parameter is not available and is set to zero. {Field 16} 

AXIS 

This 4-byte real variable is the major semi-axis for the current orbit. It is in units of 
kilometers and is denoted by a in Section 6. For T A tapes prior to September 1, 1989, this 
orbit parameter is not available and is set to zero. {Field 17} 
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ECC 

This 4-byte real variable is the eccentricity for the current orbit. It is a pure number with 
no units and is denoted by e in Section 6. For T A tapes prior to September 1, 1989, this 
orbit parameter is not available and is set to zero. {Field 18} 

ANGPER 

This 4- byte real variable is the angle of perigee for the current orbit. It is in units of de- 
grees and is denoted by T 0 in Section 6. For T A tapes prior to September 1, 1989, this or- 
bit parameter is not available and is set to zero. {Field 19} 

ISAT 

This 4-byte integer is the spacecraft identification number. For example, ISAT=8 denotes 
the F08 spacecraft, and ISAT=10 denotes the F10 spacecraft. {Field 3}. 

SPACERt61 

This 4-byte real array is filled with zeros. It is included in order to keep the size of the 
common /OUTDAT/ that same as for the first version of DECODE. 

ICOLDA(5.7} 

This 4-byte integer array contains the cold calibration counts C c for the SSM/I A-scan. 

The first dimension of the array denotes the five samples of counts that are taken for each 

channel. The second dimension of the array denotes the channel, with the order being 
19V, 19H, 22V, 37V, 37H, 85V and 85H. {Field 20} 

IHOTA15-71 

This 4-byte integer array contains the hot calibration counts C H for the SSM/I A-scan. 
The first dimension of the array denotes the five samples of counts that are taken for each 
channel. The second dimension of the array denotes the channel, with the order being 
19V, 19H, 22V, 37V, 37H, 85V, and 85H. {Field 21} 

ICOLDB15.2} 

This 4-byte integer array contains the cold calibration counts C c for the SSM/I B-scan. 

The first dimension of the array denotes the five samples of counts that are taken for each 

channel. The second dimension of the array denotes the channel, with the order being 85V 
and 85H. {Field 23} 

IHOTB15.2) 

This 4-byte integer array contains the hot calibration counts C H for the SSM/I B-scan. The 
first dimension of the array denotes the five samples of counts that are taken for each 
channel. The second dimension of the array denotes the channel, with the order being 85V 
and 85H. {Field 24} 

ALATU281 

This 4-byte real array contains the geodetic latitudes in units of degrees for the 128 cells in 
the SSM/I A-scan. Negative latitudes correspond to south latitudes, and positive values 
correspond to north latitudes. The odd elements in the array correspond to those cells for 
which observations are taken at all 7 SSM/I channels. The even elements correspond to 
those cells for which only the 85 GHz channels are sampled. The even elements are not 
compute if the User has selected I85GHZ=0. {Field 25} 
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ALONG 28) 

This 4-byte real array contains the longitudes for the 128 cells in the SSM/I A-scan. The 
longitudes are in units of degrees east of the prime meridian. The odd elements in this ar- 
ray correspond to those cells for which observations are taken at all 7 SSM/I channels. The 
even elements correspond to those cells for which only the 85 GHz channels are sampled. 
The even elements are not computed if the User has selected I85GHZ=0. {Field 26} 

BLATU281 

This 4-byte real array contains the geodetic latitudes in units of degrees for the 128 cells in 
the SSM/I B-scan. Negative latitudes correspond to south latitudes, and positive values 
correspond to north latitudes. On the B-scan, only the 85 GHz channels are sampled. This 
array is not computed if the User has selected I85GHZ=0. {Fields 25 and 27} 

BLONG28) 

This 4-byte real array contains the longitudes for the 128 cells in the SSM/I B-scan. The 
longitudes are in units of degrees east of the prime meridian. On the B-scan, only the 
85 GHz channels are sampled. This array is not computed if the User has selected 
I85GHZ=0. {Fields 26 and 27} 

TALOC5.64) 

This 4-byte real array contains the antenna temperatures (or brightness temperatures if the 
User has selected ITB=1) in units of Kelvin for the SSM/I five lower channels. The first 
dimension ich denotes the channel, with the order being 19V, 19H, 22V, 37V, and 37H for 
ich = 1 to 5. The second dimension denotes the cell position on the A-scan. The 64 cell 
positions corresponds to the odd elements in the ALAT and ALON arrays. For example, 
the latitudes for TALO(/c/r,l), TALO(/c/i,32), and TALO(/cA,64) are given by ALAT(l), 
ALAT(63), and ALAT(127), respectively. {Field 28} 

ATAHII2.128) 

This 4- byte real array contains the antenna temperatures (or brightness temperatures if the 
User has selected ITB=1) in units of Kelvin for the 85 GHz channels on the A-scan. The 
first dimension denotes the 85V and 85H channel, respectively. The second dimension 
denotes the cell position. This array is not computed if the User has selected I85GHZ=0. 
{Field 29} 

BTAHII2.128) 

This 4-byte real array contains the antenna temperatures (or brightness temperatures if the 
User has selected ITB=1) in units of Kelvin for the 85 GHz channels on the B-scan. The 
first dimension denotes the 85V and 85H channel, respectively. The second dimension 
denotes the cell position. This array is not computed if the User has selected I85GHZ=0. 
{Field 29} 

IATOILU28! 

This 4-byte integer array contains the surface-type indices for the 128 cells in the SSM/I 
A-scan. Table 5 gives the definition of the surface-type indices. These indices are com- 
puted during the FNOC processing and correspond to the FNOC latitudes and longitudes 
rather than our recomputed latitudes and longitudes that are stored in ALAT and ALON. 
The difference between the FNOC locations and our locations rarely exceeds 25 km, and 
hence IATOIL will be correct except within 25 km of boundaries. The order of this array 
is the same as ALAT and ALON. The odd elements in this array correspond to those cells 
for which observations are taken at all 7 SSM/I channels. The even elements correspond to 
those cells for which only the 85 GHz channels are sampled. The even elements are not 
computed if the User has selected I85GHZ=0. {Field 28} 
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IBTOILn281 

This 4-byte integer array contains the surface-type indices for the 128 cells in the SSM/I 
B-scan. Table 5 gives the definition of the surface-type indices. These indices are com- 
puted during the FNOC processing and correspond to the FNOC latitudes and longitudes 
rather than our recomputed latitudes and longitudes that are stored in BLAT and BLON. 
The difference between the FNOC locations and our locations rarely exceeds 25 km, and 
hence IBTOIL will be correct except within 25 km of boundaries. The order of this array 
is the same as BLAT and BLON. This array is not computed if the User has selected 
I85GHZ=0. {Field 28} 


Table 5. Definition of Surface-Type Index 


surface-type index 

definition 

0 

land 

1 

vegetal-covered land 

2 

not used 

3 

permanent sea ice 

4 

possible sea ice (note 1) 

5 

water 

6 

coast 

7 

index not available 


Note 1. Possible sea ice is a climatology of regions 
where sea ice may occur. 


Note 2. There can be a 25-km error in specifying 
surface type. 

Note 3. The surface-type indices have no seasonai 
dependence. 

Note 4. A value of 7 indicates that the FNOC surface 
index was out of bounds, and thus no surface index 
was available. 
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APPENDIX A. LIST OF ERRONEOUS DATA PERIODS 
FOR F08 SATELLITE FROM JULY 1987 THROUGH JUNE 1991 

(File BADLOC on enclosed 3.5" diskette) 
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1987 198 4.0 1987 198 5.0 
1987 217 21.5 1987 217 23.5 
1987 219 12.0 1987 219 14.0 
1987 228 14.0 1987 228 16.0 
1987 233 4.5 1987 233 6.5 
1987 236 13.5 1987 239 15.0 
1987 246 13.5 1987 246 14.5 
1987 250 11.0 1987 250 12.0 
1987 252 14.0 1987 252 16.0 
1987 253 10.5 1987 253 12.5 
1987 255 13.5 1987 255 15.0 
1987 263 16.5 1987 263 18.0 
1987 267 21.5 1987 267 23.0 
1987 269 14.0 1987 269 16.0 
1987 272 14.5 1987 272 16.0 
1987 273 15.0 1987 273 15.5 
1987 277 14.0 1987 277 15.5 
1987 280 11.5 1987 280 13.0 
1987283 1.01987283 3.0 
1987 284 21.3 1987 284 23.0 
1987 294 21.0 1987 294 22.0 
1987 304 22.0 1987 304 23.7 
1987 310 14.0 1987 310 15.0 
1987 312 5.0 1987 312 5.7 
1987 316 21.0 1987 316 23.5 
1987 317 14.2 1987 317 16.0 
1987 321 5.3 1987 321 7.8 
1987 323 5.3 1987 323 6.0 
1987 329 2.0 1987 329 3.7 

1987 329 10.1 1987 329 10.8 

1988 018 2.5 1988 018 4.5 
1988 018 9.0 1988 018 10.3 
1988 019 15.8 1988 019 17.5 
1988 022 5.3 1988 022 8.6 
1988 025 4.8 1988 025 7.2 
1988 026 10.8 1988 026 12.6 
1988 026 23.0 1988 027 1.0 
1988 030 22.1 1988 030 23.5 
1988 031 11.5 1988 031 13.3 
1988 032 22.0 1988 032 23.3 
1988 035 6.1 1988 035 6.8 
1988 067 4.5 1988 067 5.1 
1988 069 14.1 1988 069 15.8 
1988 075 23.3 1988 076 0.2 
1988 076 10.9 1988 076 13.1 
1988 084 7.5 1988 084 9.8 
1988 091 21.6 1988 092 1.1 
1988 092 21.4 1988 092 23.5 
1988 108 1.1 1988 108 2.0 
1988 110 0.7 1988 110 3.3 
1988 110 12.2 1988 110 16.0 
1988 114 13.2 1988 114 15.2 
1988 117 4.2 1988 117 7.5 
1988 118 0.7 1988 118 2.3 
1988 118 13.9 1988 118 16.0 
1988 120 5.2 1988 120 6.9 
1988 138 23.8 1988 139 1.0 
1988 140 1.2 1988 140 3.6 
1988 145 3.7 1988 145 5.6 
1988 145 15.2 1988 145 18.0 
1988 146 5.1 1988 146 7.5 
1988 147 23.5 1988 148 2.0 
1988 148 11.2 1988 148 12.5 
1988 148 14.9 1988 148 16.4 
1988 149 17.2 1988 149 21.1 
1988 149 23.0 1988 150 0.5 
1988 150 10.9 1988 150 11.8 
1988 152 16.5 1988 152 17.9 
1988 153 5.5 1988 153 9.1 
1988 154 20.6 1988 154 22.4 
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1988 

166 

21.2 

1988 

166 

22.9 

1988 

169 

2.2 

1988 

169 

4.2 

1988 

169 

10.4 

1988 

169 

11.0 

1988 

171 

15.2 

1988 

171 

15.4 

1988 

176 

14.0 

1988 

176 

15.7 

1988 

180 

19.9 

1988 

180 

21.0 

1988 

180 

23.6 

1988 

181 

0.6 

1988 

181 

1.2 

1988 

181 

1.6 

1988 

181 

4.9 

1988 

181 

6.1 

1988 

181 

9.9 

1988 

181 

11.2 

1988 

181 

20.7 

1988 

181 

21.0 

1988 

182 

3.6 

1988 

182 

4.0 

1988 

182 

5.1 

1988 

182 

5.3 

1988 

182 

12.9 

1988 

182 

13.5 

1988 

182 

17.8 

1988 

182 

19.7 

1988 

188 

11.5 

1988 

188 

13.2 

1988 

192 

14.1 

1988 

192 

16.0 

1988 

217 

20.0 

1988 

217 

21.0 

1988 

224 

13.6 

1988 

224 

15.3 

1988 

225 

9.0 

1988 

225 

11.4 

1988 

225 

12.5 

1988 

225 

13.1 

1988 

227 

3.8 

1988 

227 

5.2 

1988 

229 

16.3 

1988 

229 

16.8 

1988 

230 

12.0 

1988 

230 

13.7 

1988 

230 

16.4 

1988 

230 

19.5 

1988 

230 

21.5 

1988 

230 

22.8 

1988 

231 

15.7 

1988 

231 

16.1 

1988 

232 

1.1 

1988 

232 

3.3 

1988 

233 

17.3 

1988 

233 

21.3 

1988 

234 

12.2 

1988 

234 

12.8 

1988 

236 

10.1 

1988 

236 

11.2 

1988 

236 

11.9 

1988 

236 

12.6 

1988 

237 

22.1 

1988 

237 

22.8 

1988 

237 

23.7 

1988 

238 

0.9 

1988 

238 

2.1 

1988 

238 

2.6 

1988 

239 

21.7 

1988 

239 

23.6 

1988 

240 

0.0 

1988 

240 

0.5 

1988 

240 

2.0 

1988 

240 

4.0 

1988 

240 

9.2 

1988 

240 

10.0 

1988 

241 

7.7 

1988 

241 

8.0 

1988 

241 

9.0 

1988 

241 

12.2 

1988 

242 

4.1 

1988 

242 

6.2 

1988 

242 

7.7 

1988 

242 

8.0 

1988 

242 

21.1 

1988 

242 

22.4 

1988 

242 

22.9 

1988 

242 

23.3 

1988 

243 

3.9 

1988 

243 

4.6 

1988 

243 

9.8 

1988 

243 

10.0 

1988 

243 

11.6 

1988 

243 

12.6 

1988 

243 

16.2 

1988 

243 

19.1 

1988 

243 

21.1 

1988 

243 

21.4 

1988 

244 

0.1 

1988 

244 

1.0 

1988 

244 

2.4 

1988 

244 

5.5 

1988 

244 

9.0 

1988 

244 

9.8 

1988 

244 

15.3 

1988 

244 

16.6 

1988 

246 

6.9 

1988 

246 

8.0 

1988 

246 

11.7 

1988 

246 

12.2 

1988 

246 

14.9 

1988 

246 

16.4 

1988 

249 

10.2 

1988 

249 

12.0 

1988 

249 

12.7 

1988 

249 

13.6 

1988 

249 

17.6 

1988 

249 

20.3 

1988 

250 

4.1 

1988 

250 

7.6 

1988 

250 

8.9 

1988 

250 

9.6 

1988 

250 

12.4 

1988 

250 

13.2 

1988 

250 

15.6 

1988 

250 

15.8 

1988 

251 

15.6 

1988 

251 

16.1 

1988 

254 

9.8 

1988 

254 

15.8 

1988 

254 

17.4 

1988 

254 

18.6 

1988 

255 

3.1 

1988 

255 

4.1 

1988 

255 

14.8 

1988 

255 

15.4 

1988 

255 

18.0 

1988 

255 

21.9 
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1988 

255 

23.3 

1988 

255 

23.5 

1988 

256 

3.8 

1988 

256 

4.1 

1988 

256 

11.1 

1988 

256 

12.7 

1988 

257 

7.6 

1988 

257 

9.1 

1988 

257 

10.9 

1988 

257 

14.4 

1988 

258 

0.8 

1988 

258 

2.2 

1988 

258 

8.9 

1988 

258 

10.0 

1988 

258 

10.8 

1988 

258 

11.9 

1988 

258 

15.5 

1988 

258 

15.7 

1988 

258 

16.7 

1988 

258 

19.8 

1988 

258 

20.9 

1988 

258 

21.1 

1988 

259 

22.7 

1988 

259 

23.7 

1988 

260 

2.2 

1988 

260 

3.3 

1988 

260 

10.3 

1988 

260 

10.6 

1988 

260 

13.8 

1988 

260 

14.2 

1988 

260 

14.7 

1988 

260 

14.9 

1988 

261 

22.8 

1988 

261 

23.0 

1988 

262 

22.9 

1988 

262 

23.4 

1988 

263 

2.6 

1988 

263 

3.4 

1988 

263 

9.7 

1988 

263 

10.4 

1988 

263 

10.9 

1988 

263 

11.4 

1988 

263 

23.0 

1988 

263 

23.2 

1988 

264 

2.6 

1988 

264 

2.8 

1988 

264 

4.7 

1988 

264 

6.6 

1988 

265 

7.7 

1988 

265 

8.3 

1988 

265 

8.6 

1988 

265 

9.1 

1988 

265 

11.3 

1988 

265 

12.7 

1988 

265 

23.1 

1988 

266 

0.8 

1988 

284 

15.7 

1988 

284 

15.9 

1988 

284 

16.9 

1988 

284 

20.2 

1988 

284 

20.7 

1988 

284 

20.9 

1988 

286 

8.3 

1988 

286 

8.9 

1988 

286 

15.0 

1988 

286 

16.2 

1988 

286 

16.6 

1988 

286 

21.2 

1988 

286 

22.2 

1988 

286 

23.4 

1988 

289 

7.7 

1988 

289 

8.3 

1988 

289 

9.3 

1988 

289 

11.0 

1988 

289 

11.7 

1988 

289 

12.0 

1988 

289 

14.5 

1988 

289 

15.2 

1988 

290 

12.7 

1988 

290 

15.1 

1988 

290 

15.4 

1988 

290 

16.0 

1988 

290 

16.8 

1988 

290 

17.3 

1988 

290 

17.8 

1988 

290 

22.3 

1988 

290 

23.0 

1988 

290 

23.2 

1988 

291 

0.8 

1988 

291 

1.6 

1988 

291 

8.9 

1988 

291 

12.0 

1988 

291 

12.6 

1988 

291 

13.5 

1988 

291 

21.2 

1988 

291 

22.5 

1988 

292 

0.5 

1988 

292 

0.9 

1988 

292 

7.5 

1988 

292 

8.5 

1988 

292 

8.8 

1988 

292 

10.1 

1988 

292 

10.4 

1988 

292 

11.5 

1988 

292 

13.2 

1988 

292 

13.5 

1988 

292 

14.5 

1988 

292 

15.2 

1988 

292 

17.0 

1988 

292 

17.3 

1988 

292 

17.7 

1988 

292 

18.3 

1988 

292 

19.0 

1988 

292 

19.8 

1988 

292 

21.0 

1988 

292 

22.6 

1988 

293 

5.6 

1988 

293 

7.5 

1988 

293 

16.8 

1988 

293 

18.4 

1988 

299 

7.3 

1988 

299 

7.8 

1988 

301 

16.9 

1988 

301 

18.8 

1988 

301 

20.8 

1988 

301 

21.5 

1988 

301 

22.9 

1988 

301 

23.5 

1988 

303 

9.9 

1988 

303 

10.1 

1988 

303 

15.0 

1988 

303 

16.6 

1988 

307 

18.0 

1988 

307 

18.3 

1988 

308 

8.8 

1988 

308 

10.6 

1988 

311 

2.0 

1988 

311 

3.7 

1988 

311 

5.6 

1988 

311 

5.8 
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1988 311 23.7 1988 312 0.0 
1988 312 3.1 1988 312 5.1 
1988 312 5.9 1988 312 7.5 
1988 312 9.7 1988 312 10.0 
1988 312 11.0 1988 312 11.6 
1988 312 16.4 1988 312 17.3 
1988313 1.31988313 1.5 
1988 313 4.7 1988 313 6.2 
1988 313 6.5 1988 313 7.7 
1988 313 9.5 1988 313 9.8 
1988 314 15.5 1988 314 15.8 
1988 315 14.0 1988 315 15.8 
1988 320 17.6 1988 320 17.9 
1988 321 3.2 1988 321 4.7 
1988 321 8.9 1988 321 13.5 
1988 321 14.3 1988 321 14.8 
1988 322 9.3 1988 322 10.4 
1988 322 12.8 1988 322 14.3 
1988 325 16.9 1988 325 18.1 
1988 326 2.0 1988 326 3.2 
1988 326 4.9 1988 326 5.2 
1988 326 22.4 1988 326 22.9 
1988327 0.01988327 1.8 
1988 327 5.3 1988 327 6.1 
1988 327 11.8 1988 327 12.8 
1988 327 20.7 1988 327 21.9 
1988328 2.01988328 3.1 
1988 328 18.3 1988 328 18.9 
1988 329 17.9 1988 329 20.0 
1988 330 23.1 1988 331 1.0 
1988 331 9.1 1988 331 10.8 
1988 332 12.4 1988 332 13.8 
1988 332 14.1 1988 332 15.2 
1988 332 18.9 1988 332 19.9 
1988 334 5.6 1988 334 7.8 
1988 335 20.7 1988 335 21.3 
1988 336 5.1 1988 336 7.4 
1988 336 8.0 1988 336 9.1 
1988 336 16.5 1988 336 17.9 
1988 337 11.8 1988 337 12.8 
1988 338 7.7 1988 338 9.2 
1988 338 11.1 1988 338 11.7 
1988 338 18.1 1988 338 19.5 
1988 338 20.0 1988 338 20.8 
1988 338 21.9 1988 338 22.3 
1988339 0.41988339 1.4 
1988 339 2.7 1988 339 3.8 
1988 339 4.3 1988 339 4.9 
1988 339 14.3 1988 339 15.0 
1988 349 3.9 1988 349 4.6 
1988 349 9.4 1988 349 10.9 
1988 349 12.3 1988 349 12.5 
1988 349 12.8 1988 349 14.9 
1988 349 15.6 1988 349 16.8 
1988 349 17.0 1988 349 18.1 
1988 349 18.5 1988 349 20.7 
1988 349 21.0 1988 349 22.5 
1988 349 22.7 1988 349 23.4 
1988 350 3.9 1988 350 5.2 
1988 350 5.4 1988 350 5.7 
1988350 6.01988350 8.1 
1988 350 8.6 1988 350 9.9 
1988 350 10.2 1988 350 11.5 
1988 350 12.0 1988 350 12.4 
1988 350 16.8 1988 350 18.0 
1988 351 16.6 1988 351 19.7 
1988 351 20.0 1988 351 20.3 
1988 351 21.6 1988 351 21.8 
1988 351 23.8 1988 352 0.7 
1988352 1.61988352 2.3 
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1988 

352 

2.9 

1988 

352 

3.1 

1988 

352 

4.8 

1988 

352 

5.3 

1988 

352 

5.6 

1988 

352 

7.1 

1988 

352 

8.1 

1988 

352 

8.9 

1988 

352 

9.8 

1988 

352 

10.6 

1988 

354 

9.8 

1988 

354 

11.0 

1988 

354 

11.2 

1988 

354 

11.9 

1988 

354 

16.2 

1988 

354 

17.5 

1988 

356 

6.0 

1988 

356 

6.6 

1988 

357 

2.7 

1988 

357 

3.0 

1989 

7 

5.8 

1989 

7 

8.8 

1989 

30 

14.1 

1989 

30 

15.8 

1989 

46 

14.2 

1989 

46 

15.8 

1989 

56 

3.2 

1989 

56 

4.2 

1989 

56 

4.8 

1989 

56 

7.4 

1989 

92 

1.4 

1989 

92 

2.0 

1989 

130 

23.5 

1989 

131 

0.0 

1989 

155 

23.7 

1989 

156 

0.0 

1989 

156 

23.4 

1989 

157 

0.0 

1989 

163 

23.6 

1989 

164 

0.0 

1989 

164 

23.3 

1989 

165 

0.0 

1989 

169 

22.2 

1989 

170 

0.0 

1989 

173 

23.1 

1989 

174 

0.0 

1989 

179 

23.6 

1989 

180 

0.0 

1989 

182 

22.9 

1989 

183 

0.0 

1989 

186 

23.7 

1989 

187 

0.0 

1989 

187 

23.5 

1989 

188 

0.0 

1989 

195 

23.6 

1989 

196 

0.0 

1989 

207 

22.7 

1989 

208 

0.0 

1989 

234 

23.6 

1989 

235 

0.0 

1989 

271 

20.8 

1989 

271 

22.5 

1989 

276 

23.0 

1989 

277 

1.0 

1990 

019 

16.6 

1990 

019 

17.0 

1990 

065 

23.9 

1990 

066 

0.0 

1990 

169 

3.1 

1990 

169 

5.0 

1990 

216 

14.2 

1990 

216 

16.2 

1990 

243 

9.9 

1990 

243 

11.4 

1990 

272 

22.1 

1990 

273 

0.1 

1990 

282 

9.7 

1990 

282 

10.4 

1990 

294 

0.0 

1990 

294 

3.0 

1990 

300 

20.5 

1990 

300 

21.5 

1990 

301 

6.0 

1990 

301 

7.5 

1990 

307 

18.7 

1990 

307 

19.8 

1990 

308 

0.0 

1990 

308 

1.0 

1990 

309 

0.0 

1990 

309 

1.0 

1991 

74 

0.0 

1991 

74 

0.8 

1991 

98 

18.5 

1991 

98 

19.3 

1991 

134 

15.0 

1991 

134 

16.0 

1991 

181 

22.1 

1991 

181 

22.6 



> 
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APPENDIX B. LISTING OF SUBROUTINE DECODE, REVISION 1 

(File DECODE. FOR on the enclosed 3.5" diskette) 
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SUBROUTINE DECOOE( I85GHZ, ITS, IAOJ, IREC, LREC) 


••*•***.*••*•*•••* DESCRIPTION of arguments ************************** 
I85GHZ=0 PROCESSES ONLY THE 19, 22, AND 37 GHZ OBS. 

I85GHZ=1 PROCESSES ALL CHANNELS 

ITB=0 COMPUTES TA'S 
ITB=1 COMPUTES TB'S 

IADJ=0 DOES NO ALONGTRACK LOCATION CORRECTION 
IADJ=1 DOES AN ALONGTRACK LOCATION CORRECTION 

IREC DENOTES LOGICAL RECORD NUMBER (1 THROUGH 16) 

LREC IS THE ARRAY THAT CONTAINS THE BYTES IN THE TAPE BLOCK 
THE USER MUST SUPPLY THE TAPE READ ROUTINE FOR FILLING LREC 
BEFORE CALLING THIS ROUTINE 

************ DECLARATION OF VARIABLES, ARRAYS, AND COMMONS *********** 

CHARACTERS LREC(1784,16) 

INTEGER*2 IBUF 

COMMON /INDATA/ IBUF (1784) 

REAL*8 REV.XTIME 

INTEGER*4 ITIME, ITIMSC, IVOLT, IAGC, ICOLDA, IHOTA, ICOLDB, IHOTB, IASCTM 
INTEGER*4 IATOIL, IBTOIL, ISAT 

REAL*4 XLATSC,XLONSC,ALTSC,THT, HLTEMP,RFTEMP, FRTEHP 
REAL*4 PER I OD , ASCLOC , ANG I NC, AX I S , ECC , ANGPER , SPACER 
REAL*4 ALAT ,ALON,BLAT,BLON,TALO,ATAHI ,BTAHI 

COMMON/OUTDAT/ REV, XT I ME , I T I ME , I T I MSC , XLATSC, XLONSC, ALTSC , THT , 

1 HLTEMP(3), 1V0LT(2),RFTEMP, FRTEMP, IAGC(6), 

2 IASCTM, PERIOO, ASCLOC, ANGINC, AXIS, ECC, ANGPER, ISAT, SPACER(6), 

2 ICOLDA(5,7), IHOTA(5,7), ICOLDB(5,2), IHOTB(5,2), 

3 ALAT(128),ALON(128),BLAT(128),BLON(128), 

4 TALO(5,64),ATAHI(2, 128),BTAHI (2, 128), IATOILC128), IBT0ILC128) 

********************** data initialization *********************** 


DATA N1,N2, 143/256,65536, 16777216/ 

DATA RAD/0.017453293/ 

************************ gEQjfj EXECUTION **************************** 


TRANSFER TAPE BYTES IN IREC LOGICAL RECORD TO ARRAY IBUF 

DO 100 1=1,1784 

IBUFCI )=ICHAR(LREC( I , IREC)) 

100 CONTINUE 
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FIND SCAN TIME 

ITIME= N3*IBUF( 1 )+N2*18UF( 2)+N1*lBUF< 3)+IBUF< 4) 
IFRCTM=N3*IBUF(17)+N2*IBUF(18)+N1*I8UF(19)+IBUF(20) 
IF(IFRCTM.EQ.O) XTIME=IT1ME 

IF(IFRCTM.NE.O) XTIME=ITIME+1 .D-4*(I FRCTM-10000) 

FIND SPACECRAFT TIME AND ORBIT NUMBER 
IOLD-1 FOR ALL SCANS BEFORE 1989 

IFdTIME.LT. 63163966) THEN 
IOLD-1 

ITIMSC= N3*I8UF( 9)+N2*IBUF<10)+N1*IBUF<11)+IBUF<12) 

ELSE 

IOLD s O 

ITIMSC-ITIME 
END IF 

IF(ITIME. GE.63163966.AND.ITIME.LT. 84156110) THEN 
REV=1.0-4*(N3*IBUF( 9)+N2*lBUF(10)+N1*IBUF(1 1 )+IBUF(12)) 

ELSE 

REV=1.D-4*(N3*IBUF( 5)+N2*IBUF< 6)+N1*IBUF( 7)+IBUF( 8)) 

ENDIF 

FIND SPACECRAFT LATITUDE, LONGITUDE, AND ALTITUDE 

XLATSC=1.D-6*(N3*IBUF(13)+N2*IBUF<14)+N1*IBUF(15)+IBUF(16))-90. 
XLONSC=1.D-6*(N3*IBUF(21)+N2*IBUF(22)+Nl*IBUF(23)+IBUF(24)) 
ALTSC= 1 .D-3*(N3*IBUF(25)+N2*IBUF(26)+N1*1BUF(27)+IBUF(28)) 

DETERMINE SATELLITE NUMBER AND INCIDENCE ANGLE 
FACTOR OF .7040147 = SIN(44.75 DEG) 

FACTOR OF .7116583 = SIN<45.37 DEG) 

IFdTIME.GT. 144554200) THEN 

ITERM=N3*IBUF( 9)+N2*IBUF(10)+N1*IBUF(11)+IBUF(12) 
THT=0.001*INT(ITERM/1000) 

I SAT= I TERM-1 000* I NT (I TERM/ 1 000 ) 

ELSE 

EARTHR=RADIUS OF CURVATURE OF EARTH'S SURFACE 
EARTHR=6345 . 7+55 . 0*S I N ( RAD*XLATSC)**2 
D I FREV=DABS< REV-300- (XT I ME- 1 6530609 )/61 18.) 

I F(D IFREV.LT. 100.) THEN 
ISAT=8 

THT=ASIN(0.7040147*(EARTHR+ALTSC)/EARTHR)/RAD 

ELSE 

ISAT=10 

THT=ASIN(0.7116583*(EARTHR+ALTSC)/EARTHR)/RAD 

ENDIF 

ENDIF 


FIND RADIOMETER CALIBRATION DATA 

• HLTEMP(3)=0.01*(N1*IBUF(29)+IBUF(30)) 
HLTEMP(2)=0.01*(N1*IBUF(31)+IBUF(32)) 
HLTEMP(1)=0.01*(N1*IBUF(33)+IBUF(34)) 
RFTEMP= 0.01*(N1*IBUF(39)+IBUF(40)) 
FRTEMP= 0.01*(N1*IBUF(41)+IBUF(42)) 
IVOLT(2)= N1*IBUF(35)+I8UF<36) 

IVOLT(1)= N1*IBUF(37)+IBUF<38) 

IAGC(3)= N1*IBUF(43)+IBUF(44) 

IAGC(2)= N1*I8UF(45)+IBUF(46) 

IAGC(1 )= N1*IBUF(47)+IBUF(48) 
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N=-2 

DO 200 ICH=1,7 
IF(ICH.NE.7) SPACER(ICH)=0. 

DO 200 IP=1 ,5 
N=N+2 

I COLD A( IP, ICH)=N1*IBUF(N + 77)+IBUF(N + 78) 

IH0TA( IP, ICH)= N1*IBUF(N +147)+IBUF(N +148) 

200 CONTINUE 
C 

1AGC(6)=N1*IBUF(217)+IBUF(218) 

IAGC(5)=N1*IBUF(219)+I8UF(220) 

IAGC(4)=N1*IBUF(221 )+IBUF(222) 

C 

N=-2 

DO 300 ICH=1 ,2 
DO 300 IP=1 ,5 
N=N+2 

IC0LDB( IP, ICH)=N1*IBUF(N +223)+IBUF(N +224) 

IH0TB( IP, ICH)= N1*IBUF(N +243)+IBUF(N +244) 

300 CONTINUE 

FIND ORBIT PARAMETERS 

I F( ITIME.LT. 841 561 10) THEN 
IASCTM=0 
PERIOO=0. 

ASCLOC=0. 

ANGINC=0. 

AXIS=0. 

ECC=0. 

ANGPER-0. 

ELSE 

IASCTM= N3*IBUF(49)+N2*IBUF(50)+N1*1BUF(51 )+IBUF(52) 

PERIOD 3 1 .D-4*(N3*IBUF(53)+N2*IBUF(54)+N1*IBUF(55)+I8UF(56)) 
ASCLOC= 1 .D-3*(N3*IBUF(57)+N2*IBUF(58)+N1*IBUF(59)+IBUF(60)) 
ANGINC= 1 .D-6*(N3*IBUF(61 )+N2*IBUF(62)+N1*IBUF(63)+IBUF(64)) 
ANG I NC=1 80 . -ANG I NC 

AXIS 3 1 .D-4*(N3*IBUF(65)+N2*IBUF(66)+N1*IBUF(67)+IBUF(68)) 
ECC= 1.D-10*(N3*5BUF(69)+N2*!BUF(70)+N1*IBUF(71)+IBUF(72)) 
ANGPER= 1.D-5*(N3*IBUF(73)+N2*IBUF(74)+N1*IBUF(75)+IBUF(76)) 
END IF 

DETERMINE IF ALONGTRACK AND YAW ADJUSTMENT IS NEEDED 
TRKADJ=-1 .2 CORRESPONDS TO A 15 KM ALONGTRACK ADJUSTMENT 
YAW 3 . 3125 CORRESPONDS TO A 0.5 DEG. YAW ROTATION 

TRKADJ 3 0. 

IF( IOLD.EQ. 1 .AND. IADJ.EQ.1 ) TRKADJ=-1.2 
YAW=0 . 

IF(ISAT.EQ.8) YAW 3 . 3 125 

FIND CELL LATITUDES AND LONGITUDES 

J85GHZ 3 I85GHZ 
IF(TRKADJ.NE.O.) J85GHZ=1 
CALL FDLTLN( J85GHZ , I OLD , YAW) 

IF(TRKADJ.NE.O.) CALL ADJLOC(TRKADJ) 

FIND EITHER TA'S OR TB'S 

CALL FDTA( I85GHZ) 

I F( ITB.EQ. 1 .AND. ISAT.EQ. 8) CALL FDTB08(I85GHZ) 

IF( ITB.EQ. 1 .AND. I SAT. NE. 8) CALL FDTB00(I85GHZ) 

RETURN 
END 
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SUBROUTINE FOLTLN( J85GHZ, I OLD, YAW) 

C 

C THIS SUBROUTINE FINDS THE LATS AND LONS FOR THE SSMI CELLS 
C J85GHZ=0 DOES NOT DO 85 GHZ LAT/LON, J85GHZ=1 DOES 85 GHZ LAT LON 
C 

INTEGER*4 INDEX(19), JNDEX(3,109) 

C 

C SPECIFY COMMON /INDATA/ 

C 

INTEGER*2 IBUF 

COMMON /INDATA/ IBUF(1784) 

C 

C SPECIFY COMMON /OUTDAT/ 

C 

REAL*8 REV.XTIME 

INTEGER*4 ITIME, ITIMSC, IVOLT, IAGC, ICOLDA, IHOTA, ICOLDB, IHOTB, IASCTM 
INTEGER*4 IATOIL, IBTOIL, I SAT 

REAL*4 XLATSC,XLONSC,ALTSC,THT,HLTEMP,RFTEMP, FRTEMP 

REAL*4 PERIOD, ASCLOC, ANGI NC , AX I S , ECC , ANGPER , SPACER 

REAL*4 ALAT , ALON , BLAT , BLON , TALO, ATAH I , BTAH I 

COMMON/OUTDAT/ REV, XTIME, ITIME, ITIMSC, XLATSC,XLONSC,ALTSC,THT, 

1 HLTEMP(3), I VOLT(2),RFTEMP, FRTEMP, IAGC(6), 

2 IASCTM, PERIOD, ASCLOC, ANGINC, AXIS, ECC, ANGPER, ISAT,SPACER(6), 

2 ICOLDA(5,7), IHOTA(5,7), ICOLDB(5,2), IHOTB(5,2), 

3 ALAT(128),AL0N(128),BLAT(128),BL0N(128), 

4 TALO(5,64),ATAHI (2, 128),BTAHI(2, 128), IATOILC128), IBT0ILO28) 

C 

C DATA INITIALIZATION 
C 

DATA N 1/256/ 

DATA RAD/0.017453293/ 

DATA INDEX/1,9,17,25,33,41,49,57,65,73,81,89,97,105,113,121,123, 

1 127,128/ 

DATA JNDEX/ 

1 5, 1, 9, 13, 9, 17, 21, 17, 25, 29, 25, 33, 37, 33, 41, 

1 45, 41, 49, 53, 49, 57, 61, 57, 65, 69, 65, 73, 77, 73, 81, 

1 85, 81, 89, 93, 89, 97,101, 97,105,109,105,113,117,113,121, 

1 3, 1, 5, 7, 5, 9, 11, 9, 13, 15, 13, 17, 19, 17, 21, 

1 23, 21, 25, 27, 25, 29, 31, 29, 33, 35, 33, 37, 39, 37, 41, 

1 43, 41, 45, 47, 45, 49, 51, 49, 53, 55, 53, 57, 59, 57, 61, 

1 63, 61, 65, 67, 65, 69, 71, 69, 73, 75, 73. 77. 79, 77. 81, 

1 83, 81, 85, 87, 85, 89, 91, 89, 93, 95. 93, 97, 99, 97,101, 

1 103,101,105,107,105,109,111,109,113,115,113,117,119,117,121, 

1 125,123,127, 2, 1. 3, 4. 3, 5, 6, 5, 7. 8, 7. 9, 

1 10, 9, 11, 12, 11, 13, 14, 13, 15. 16, 15, 17, 18, 17, 19, 

1 20, 19, 21, 22, 21, 23, 24, 23. 25, 26, 25, 27, 28, 27, 29, 

1 30. 29, 31, 32, 31, 33, 34, 33, 35, 36, 35, 37, 38, 37, 39, 

1 40, 39, 41, 42, 41, 43. 44, 43, 45, 46, 45, 47, 48, 47, 49, 

1 50. 49, 51, 52, 51, 53, 54, 53, 55, 56, 55, 57, 58, 57, 59, 

1 60, 59. 61, 62. 61, 63. 64, 63, 65, 66. 65, 67, 68, 67, 69. 

1 70, 69, 71, 72, 71, 73, 74, 73, 75, 76, 75, 77, 78, 77, 79. 

1 80, 79, 81, 82, 81, 83, 84, 83, 85, 86, 85, 87, 88, 87, 89, 

1 90, 89, 91. 92, 91, 93, 94, 93, 95, 96, 95, 97, 98, 97, 99, 

1 100, 99,101,102,101,103,104,103.105,106,105,107,108,107.109, 

1 110,109,111,112,111,113,114,113,115.116,115,117,118,117,119, 

1 120,119,121,122,121,123,124,123,125,126,125,127/ 

C 

C BEGIN EXECUTION 

C 

C SET TABLE LAT/LON FOR A- SCAN 

C 

N=-2 

DO 100 JCEL=1,19 
N=N+2 

ICEL=INDEX( JCEL) 

ALAT ( ICEL)=0.01*(N1*IBUF(N +263)+IBUF(N +264)-9000> 
ALON(ICEL)=0.01*(N1*1BUF(N +301)+IBUF(N +302)) 
IF(ALON(ICEL).GE.360.) AL0N(rCEL)=AL0N(ICEL)-360. 

100 CONTINUE 
C 
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c 

c 


SET MID-POINTS FOR A-SCAN 


NCEL=46 

I F( J85GHZ.EQ. 1 ) NCEl-109 
DO 200 JCEL S 1 ,NCEL 
ICEL=JNDEX( 1 , JCEL) 

I1=JNDEX(2, JCEL) 

I2=JNDEX(3, JCEL) 

D I FLAT =ALAT < 1 2)-ALAT (II) 

AVGLAT=0 . 5*<ALAT( I 1 )+ALAT< 12) ) 

D I FLON=ALON( I2)-ALON( 1 1 ) 

IF(DIFLON.GT. 180.) D I F LON=D I F LON-360 . 

I F (0 1 FLON . LT . - 180 . ) D I F LON=D I F L0N+360 . 

AVGLON=ALON ( II )+0 . 5*0 1 FLON 
XS0= ( 2 . *RAD*AVGLAT )**2 

X FAC= 1 . -0 . 1 66271 42*XSQ+0 . 00807934*XSQ*XSQ- 0.000151 880*XSQ*XSQ*XSQ 
ALAT ( I CEL )=AVGLAT*( 1 . +0 . 1 25*(RAD*D I FLON )**2*XFAC ) 

X=RAD* ( 90 . -ABS ( AVGLAT ) ) 

TANLAT=1./<X+X*X*X/3.) 

1 F(AVGLAT.LT.O. ) TANLAT=-TANLAT 
ALON(ICEL)=AVGLON-0.2500*RAD*DI FLAT*DI FLON*TANLAT 
1 F(ALON(ICEL) .LT . 0.) ALON(ICEL)=ALON(ICEL)+360. 

I F ( ALON ( I CE L ) . GE . 360 . ) ALON ( I CE L ) =ALON ( I CEL ) -360 . 

200 CONTINUE 

DO YAW CORRECTION 

I F(J85GHZ.EQ.O) THEN 
IF(YAU.EQ.O.) RETURN 

DO 260 ICEL=1, 127,2 
IF(ICEL.E0.127) GO TO 250 
D I FLAT=ALAT ( I CEL+2)-ALAT (I CEL) 

DIFLON=ALON(ICEL+2)-ALON(ICEL) 

I F(D IFLON.LT. -180. ) DI FLON=DI FLON+360. 

IF(DIFLON.GT. 180.) 0IFL0N=DIFL0N-360. 

250 CONTINUE 

ALAT ( I CEL )=ALAT< I CEL )+YAW*D I FLAT 
ALON (I CEL )=ALON( I CEL )+YAW*D I FLON 
I F ( ALON (I CEL ) . LT . 0.) ALON<ICEL)=ALON(ICEL)+360. 

I F ( ALON (I CE L ) . GE . 360 . ) ALON dCEL)=ALONd CEL) -360. 

260 CONTINUE 
RETURN 
ENDIF 

SET TABLE LAT/LON FOR B-SCAN 
N=-2 

DO 300 JCEL=1 , 19 
N=N*2 

ICEL=INDEX( JCEL) 

IDEL=N1*IBUF(N +339)+IBUF(N +340) 

IFdDEL.GT. 32767) IDEL=IDEL-65536 
LATDEL=(IDEL+30000)/ 1000-30 
LONDEL= I DEL+291 00- 1 000*( LATDEL+30) 

BLAT ( I CEL )=ALAT (I CEL )+0 . 01 *LATDEL 
BLON(ICEL)=ALON(ICEL)+0.01*LONDEL 
I F(BLONdCEL) .LT. 0.) BLON(ICEL)=BLON(ICEL)+360. 

I F(BLONdCEL) .GE.360. ) BLON(ICEL)-BLON(ICEL)-360. 

300 CONTINUE 
C 
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SET MID-POINTS FOR B-SCAN 

DO 400 JCEL S 1 ,109 
ICEL=JNDEX( 1 , JCEL) 

I 1 =JNDEX(2, JCEL ) 

I2=JNDEX(3, JCEL) 

D1 FLAT=BLAT ( I 2)-BLAT (II) 

AVGLAT=0.5*(BLAT ( I 1 )+BLAT( 12)) 

DI FLON=BlON( I2)-BLON( 1 1 ) 

IF(DIF10N.GT. 180.) DIFLON=D1FLON-360. 

I F<D I FLON . LT . - 180 . ) D I FLON=D I FLON+360 . 

AVGLON=BLON <11 >+0.5*01 FLON 
XSQ= ( 2 . *RAD*AVGLAT )**2 

X F AC= 1 . -0 . 1 66271 42*XSQ+0 . 00807934*XSQ*XSQ- 0.000151 880*XSQ*XSQ*XSQ 
BLAT ( ICEL)=AVGLAT*( 1 .+0. 125*(RAD*DI FLON)**2*XFAC) 

X=RAD* ( 90 . - ABS( AVGL AT ) ) 

TANLAT=1 ./(X+X*X*X/3. ) 

IF(AVGLAT.LT.O-) TANLAT=-TANLAT 

BLON( I CEL )=AVGLON-0 . 2500*RAD*D I FLAT*D I FLON+TANLAT 

I F(BLON( ICED.LT. 0.) BLON( ICEL)=BLON( ICED+360. 

I F ( BLON ( I CE L ) . GE . 360 . ) BLON(ICEL)=BLON(ICEL>-360. 

400 CONTINUE 

CORRECT CELL 128 FOR DATA BEFORE 1989 

IF(IOLD.EQ.I) THEN 
D I FLAT=ALAT ( 1 27)-ALAT (126) 

D I F LON = ALON ( 1 27)-AL0N (126) 

I F (D I FLON . LT . -1 80 . ) D I FLON=D I FLON+360 . 

IF(DIFLON.GT. 180.) DIFLON=DIFLON-360. 

ALAT ( 1 28)=ALAT ( 1 27)+D I FLAT 
ALON ( 1 28) =ALON ( 1 27)+D I FLON 
IF(ALON(128).LT. 0.) ALON(128)=ALON(128)+360. 

I F ( ALON ( 1 28) . GE . 360 . ) ALON ( 1 28) -ALON ( 1 28 >-360 . 

D I FLAT=BLAT ( 1 27) -BLAT ( 1 26) 

D I F LON=BLON ( 1 27 ) -BLON ( 1 26 ) 

IF(DI FLON. LT. -180.) DIFL0N=D! FLON+360. 

IF(DI FLON.GT. 180.) DI FLON=DI FL0N-360. 

BLAT ( 1 28)=BLAT ( 127)+D I FLAT 
BLON(128)=BLON(127)+DI FLON 
I F(BL0N(128).LT . 0.) BLON(128)=BLON(128)+360. 

IF(BLON(128) .GE.360. ) BLON(128)=BLON(128)-360. 

END IF 

DO YAW CORRECTION 

IF(YAU.EQ.O. ) RETURN 
DO 600 ICEL=1 ,128 
I F(ICEL.EQ. 128) GO TO 580 
DI FLAT1=ALAT(ICEL+1 )-ALAT(ICEL) 

D I FLON1=ALON( ICEL+1 )-ALON( ICEL) 

I F(DI FLON1.LT. -180.) DI FLON1=DI FLON1+360. 

I F (D I FLON1 .GT. 180.) DIFLON1=DIFLON1-360. 

DIFLAT2=BLAT( ICEL+1 )-BLAT( ICEL) 

DI FLON2=BLON( ICEL+1 )-BLON( ICEL) 

I F ( D I FLON2 . LT . - 1 80 . ) D I FLON2=D I F L0N2+360 . 

IF(DIFLON2.GT. 180.) DIFLON2=DIFLON2-360. 

580 CONTINUE 

ALAT ( I CE L ) =ALAT ( I CEL )+2 . * YAW+D I F LAT 1 
ALON ( I CEL ) =ALON ( I CE L )+2 . * YAW+D I F LON 1 
I F(ALON( ICED.LT. 0.) ALON(ICEL)=ALON( ICED+360. 

I F(ALON( ICEL). GE.360.) ALON(ICEL)=ALON(ICEL)-360. 

BLAT ( I CEL)-BLAT( I CEL )+2 . *YAW*D I FLAT2 
BLON ( I CEL)=BLON( I CEL)+2 . *YAW*D I FLON2 
IF(BLON(ICEL).LT. 0.) BLON(ICEL)=BLON(ICEL)+360. 

I F(BLON( ICEL). GE.360.) BL0N(ICEL)=BL0N(ICEL)-360. 

600 CONTINUE 
RETURN 
END 
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c 

SUBROUTINE ADJLOC(TRKADJ) 

CORRECT FOR ALONGTRACK ERROR 
SPECIFY COMMON /OUTDAT/ 

REAL*8 REV.XTIME 

I NTEGER*4 I T I ME , I T I MSC , I VOLT , I AGC, I COLDA , I HOTA , I COLDB , IHOTB, IASCTM 
INTEGERS IATOIL, IBTOIl, ISAT 

REAL*4 XLATSC,XLONSC,ALTSC,THT,HLTEMP,RFTEMP, FRTEMP 
REAL*4 PERIOO, ASCLOC.ANGINC, AXIS, ECC.ANGPER, SPACER 
RE AL*4 ALAT , ALON , BLAT , BLON , TALO , AT AH I , BT AH I 
COMMON/OUTDAT/ REV, XT I ME , I T IME , I T IMSC , XLATSC , XLONSC , ALTSC , THT , 

1 HLTEMP(3), I VOLT(2),RFTEMP, FRTEMP, IAGC(6), 

2 IASCTM, PERIOO, ASCLOC.ANGINC, AXIS, ECC.ANGPER, t$AT,SPACER<6), 

2 tCOLDA(S,7),IHOTA(5,7),ICOLDB(5,2),IHOTB(5,2), 

3 ALAT(128), ALON (128), BLAT (128), BLON (128), 

4 TALO(5,64),ATAHI(2,128),BTAHI(2,128), IAT01LO28), IBTOlL(128) 

BEGIN EXECUTION 

DO 100 1CEL»1 ,128 
OIFLAT=BLAT(ICEL)-ALAT(ICEL) 

ALAT ( I CEL )=ALAT ( I CEL J+TRKAO J*D I FLAT 
BLAT ( I CEL )=BLAT ( I CEL HTRKAD J*D 1 FLAT 
C 

DIFLON=BLON( ICEL)-ALON(ICEL) 

I F(D IFLON.LT. -180.) DIFLON=DIFLON+360. 

IF(DIFLON.GT. 180.) DIFLON=DIFLON-360. 

C 

ALON ( I CEL )=ALON ( I CEL )+TRKAD J*DI FLON 
I F(ALON( ICED.LT . 0.) ALON(ICEL)=ALON<ICEL)+360. 

I F(ALON( ICEL).GE.360. ) ALON<ICEL)=ALON( ICEL)-360. 

C 

BLON(ICEL) =BLON ( I CEL )+TRKAD J*D I FLON 
I F(BLON(ICEL).LT . 0.) BLON(ICEL)=BLON(ICEL)+360. 

I F(BLON( ICEL) .GE.360. ) BLON(ICEL)=BLON( ICED-360. 

100 CONTINUE 
RETURN 
END 
C 
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SUBROUTINE F0TA(I8SGHZ) 

THIS SUBROUTINE FINDS THE ANTENNA TEMPERATURES AND SURFACE TYPES 

INTEGERS I TAL0(5 ) , I TAH I (8) 

SPECIFY COMMON /INDATA/ 

INTEGER*2 IBUF 

COMMON /INDATA/ IBUF(1784) 

SPECIFY COMMON /OUTDAT/ 

REAL*8 REV.XTIME 

INTEGERS ITIME, ITIMSC, I VOLT, I AGC , I COLD A , IHOTA, ICOLDB, I HOTB, IASCTM 
INTEGER** IATOIL, IBTOIl, ISAT 

REAL** XLATSC,XLONSC,ALTSC,THT,HLTEMP,RFTEMP,FRTEMP 

REAL** PER 1 00 , ASCLOC , ANGI NC , AX I S , ECC , ANGPER , SPACER 

REAL** ALAT , ALON , BLAT , 8L0N , TALO , AT AH I , BTAH I 

COMMON/OUTDAT/ REV.XTIME, ITIME, ITIMSC, XLATSC,XLONSC,ALTSC,THT, 

1 HLTEMP(3),IVOLT(2),RFTEMP,FRTEMP,IAGC(6), 

2 IASCTM, PERIOD, ASCLOC, ANGINC, AXIS, ECC, ANGPER, ISAT, SPACER<6), 

2 ICOLDA(5,7),IHOTA(5,7),ICOLDB(5,2),IHOTB<5,2), 

3 ALAT ( 128) , ALON( 128) , BLAT (128), BLON (128), 

* TALO(5,6*),ATAHI(2,128),BTAHI(2,128),IATOIL(128),IBTOIL<128) 

DATA N1 ,N2/256, 65536/ 

BEGIN EXECUTION 


N=-10 

M=-12 


DO 100 ICEL=1,6* 

JCEL=2*ICEL-1 

FIND THE TA'S FOR THE 3 LOWER FREQUENCIES 
N=N+10 


IU0RK4=N2*1BUF<N ♦377)+N1*IBUF<N+378)+IBUF<N+379) 
I T AV= I NT ( I WORK*/*096 ) 

ITALO(1 )=ITAV 
ITALO(2)=I WORK4-4096* I TAV 

IWORK*=N2*IBUF(N +380)+N1*IBUF(N+381 )+IBUF(N+382) 
ITAV=INT( IUORK*/*096) 

ITALO(*)=ITAV 
I TALO( 5 )= I WORK*-*096* I TAV 

IWORK*=N2*IBUF(N +383)+N1*IBUF(N+38*)+IBUF(N+385) 
I TAV- I NT ( IWORK*/*096) 

ITALO(3)=ITAV 

SCALE THE TA'S 

DO 20 ICH-1,5 

IF(ITALO(ICH).LE.3800) THEN 
TALO( I CH , I CEL )=0 . 1 * 1 TALO( I CH ) 

ELSE 

TALOdCH, ICEL)=ITALO< ICH)-3*20 
END IF 

20 CONTINUE 

FIND THE TOIL FLAGS 


C 


I RES= I UORK4-4096* I TAV 
ITOIL1-INT(IRES/512) 
IATOI L( JCEL)=IT0IL1 


65 



non n n n o n non 


IF(I85GHZ.EQ.O) GO TO 100 
FIND 85 GHZ TOILS'S AND TA'S 
M=M+12 

I RES= I RES- 1 TOIL 1*5 12 
ITOIL2=INT( IRES/64) 

IRES=IRES-IT0IL2*64 
ITOI L3=INT (IRES/8) 

ITOIL4=IRES-ITOIL3*8 
IBT0IL(JCEL)=IT0IL2 
IAT0IL( JCEL+1 )=ITOIL3 
IBTOILCJCEL+1 )=ITOIL4 

IWORK4=N2*IBUF(M+1017)+N1*1BUF(M+1018)+IBUF(M+1019) 
I TAV= I NT ( I WORK4/4096 ) 

ITAHI (1 )=1TAV 
I T AH I ( 3 ) = I WORK4-4096* I TAV 

IWORK4=N2*IBUF(M+1020)+N1*IBUF(M+1021 )+IBUF(M+1022) 
ITAV=INT( IWORK4/4096) 

ITAHI (2)=ITAV 

ITAHI <4)=I WORK4-4096* I TAV 

IWORK4=N2*IBUF(M+1023)+N1*IBUF(M+1024)+IBUF(M+1025) 
I T AV= I NT ( I WORK4/4096 ) 

ITAHI<5)=ITAV 

I T AH I ( 7)= I WORK4-4096* I TAV " 

IWORK4=N2*IBUF(M+1026)+N1*IBUF(M+1027)+IBUF(M+1028) 
ITAV=INT ( IWORK4/4096) 

ITAHI <6)=ITAV 

ITAHI (8)=I UORK4-4096* I TAV 

SCALE THE TA'S 


K=0 

DO 40 KCEL S JCEL, JCEL*1 

DO 40 ICH=1 ,2 

K=K+1 

1F(ITAHI(K).LE.3800) then 
ATAHK ICH,KCEL)=0. 1 ‘ITAHI (K) 
ELSE 

ATAHI(ICH,KCEL)=ITAHI(K)-3420 

ENOIF 

C 

K=K+1 

C 

IF(ITAHI(K).LE.3800) THEN 

BTAHI(1CH,KCEL) S 0.1*ITAHI(K) 

ELSE 

8TAHI(ICH,KCEL) S ITAHI (K)-3420 
ENDIF 
C 

40 CONTINUE 
C 

100 CONTINUE 
RETURN 
END 
C 
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SUBROUTINE FDTB08( I85GHZ) 

THIS SUBROUTINE CONVERTS ANTENNA TEMPS. TO BRIGHTNESS TEMPS. 

F08 ALONG-SCAN AND OFFSET CORRECTIONS 

REAL*4 DELTA(4),CHI(2,4),AVGTA(5),ABIAS(5) 

REAL*4 SBIAS1(64),SBIAS2(64),SBIAS3(64),SBIAS4(64),SBIAS5(64) 
REAL*4 AW(4),AHV(4),AOV(4),AHH(4).AVH(4),AOH<4) 

SPECIFY COMMON /OUTDAT/ 

REAL*8 REV, XT I ME 

INTEGERS ITIME, ITIMSC, IVOLT, IAGC, ICOLDA, IHOTA, ICOLDB, IHOTB, IASCTM 
INTEGERS IATOIL, IBTOIL, ISAT 

REAL*4 XLATSC,XLONSC,ALTSC,THT,HLTEHP,RFTEMP, FRTEMP 

REAL*4 PERIOD, ASCLOC,ANGINC, AXIS, ECC.ANGPER, SPACER 

REAL*4 ALAT , ALON , BLAT , BLON , TALO, ATAH I , BTAH I 

COMMON/OUTDAT/ REV, XT I ME, ITIME, ITIMSC,XLATSC,XLONSC,ALTSC,THT, 

1 HLTEMP(3), IVOLT(2),RFTEMP,FRTEMP, IAGC(6), 

2 I ASCTM, PER IOO, ASCLOC,ANGINC, AXIS, ECC.ANGPER, I SAT, SPACER (6), 

2 ICOLDA(5,7), IHOTA(5,7), ICOLDB(5,2), IHOTB<5,2), 

3 ALAT ( 1 28), ALON ( 1 28), BLAT ( 128), BLON( 128), 

4 TAL0(5,64),ATAHI(2,128),BTAHI(2,128), IAT0ILO28), IBTOIL(128) 

DATA INITIALIZATION 
DATA I START/ 1/ 

DATA DELTA/0.03199,0.02685,0.01434,0.01186/ 

DATA CHI/. 00379, .00525, .00983,0.0, .02136, .02664, .01387, .01967/ 

DATA AVGTA/190. 93, 130. 14,215.42,21 1 .39,158. 16/ 

********************************************************** 

DATA ABIAS/5*0./ 

IN ORDER TO APPLY WENTZ [1992] OFFSET TO TB'S, 

ACTIVATE NEXT STATEMENT AND COMMENT OUT PREVIOUS STATEMENT 
DATA ABIAS/2.0,3.5, 1 .32,-1 .57,-0.22/ 

**************** ********** ************************ ******** 


DATA SB IASI/ 





-.03, 

• 01, 

• 04, 

.06, 

+ 

.09, 

.10, 

• 12, 

.12, 

• 13. 

.13, 

• 13, 

• 13. 

•13, 

.14, 

+ 

,14, 

.14, 

.14. 

.14, 

.15. 

.15, 

.14, 

.14. 

•13, 

.12, 

+ 

■12. 

.11, 

.11, 

.11, 

.11, 

.11. 

.12. 

• 11, 

■10, 

.10, 

+ 

.09, 

.07, 

.06, 

.05, 

• 05. 

•04, 

• 03. 

.03. 

•02, 

.01, 

+ 

-.01. 

-.02, 

-.04, 

-.05, 

-.06, 

-.08, 

-.09. 

-.09, 

-.11, 

-.12. 

+ 

-.14, 

-.18, 

-.21, 

-.26, 

-.30. 

-.36. 

-.43, 

-.52, 

-.63, 

-.75/ 

DATA SBIAS2/ 





-.03, 

-.01, 

• 00, 

.01, 


.01, 

.01, 

• 01. 

• 01, 

•01, 

• 01, 

.01, 

• 02. 

• 03, 

•04, 

+ 

.04, 

•05, 

.05, 

• 05, 

•05, 

.05, 

• 04, 

.04, 

■ 04, 

•04, 

+ 

.04, 

• 04, 

• 04, 

• 04, 

.04, 

•05, 

.06, 

• 07, 

• 07, 

.08, 

+ 

.07, 

• 07, 

• 06, 

• 07, 

.07, 

•06, 

• 06, 

.06, 

.06, 

.06, 

+ 

.05, 

.04, 

• 04, 

.04, 

.04, 

• 02, 

.00, 

-.01, 

-.02, 

-.03, 

+ 

-.04, 

-.06, 

-.09, 

-.12, 

-.16, 

-.20, 

-.26, 

-.32, 

-.39, 

-.48/ 

DATA SBIAS3/ 





-.05, 

.00, 

• 04, 

.07. 

+ 

.08, 

• 09, 

■ 10. 

.11. 

.11, 

.12, 

.12, 

.13, 

■ 14, 

• 14. 

+ 

.15, 

.15, 

.17, 

.17, 

.16, 

.16, 

• 15, 

• 14. 

.14. 

• 14. 

+ 

.13, 

-13. 

.13, 

.13. 

.13. 

-13. 

.14, 

•13, 

• 14. 

•13, 

+ 

• 12. 

.11, 

.10, 

• 09, 

.09. 

•08, 

• 07, 

• 05, 

•04, 

.03, 

+ 

• 02, 

.00, 

-.01. 

-.02, 

-.04, 

-.05, 

-.07, 

-.10, 

-.12. 

-.14, 

+ 

-.17, 

-.20, 

-.25, 

-.31, 

-.36, 

-.42. 

-.50, 

-.63, 

-.79, 

-.97/ 

DATA S8IAS4/ 





-.06, 

-.03, 

o 

o 

• 01. 

+ 

.03, 

.04, 

• 05, 

.06. 

.06, 

.07. 

.09, 

.10, 

.11. 

.13, 

+ 

.14, 

• 15, 

.15, 

.16, 

.16, 

.16, 

.16, 

■ 16. 

.16, 

.15, 

+ 

.14, 

.14, 

.14, 

.13, 

.13, 

• 13. 

• 14, 

• 13, 

•13, 

■ 13, 

+ 

.12, 

.11, 

• 10, 

•10, 

.09, 

.08. 

.08, 

.06, 

.06, 

•05, 

+ 

.04, 

.03, 

• 02, 

.01, 

.00, 

-.02, 

-.05, 

-.07, 

-.09, 

-.12, 

+ 

-.16, 

-.18, 

-.22, 

-.27, 

-.33. 

-.42. 

-.50, 

-.60. 

-.74, 

-.90/ 
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DATA SBIAS5/ 





-.23, 

-.20, 

-.17, 

-.15, 

+ 

-.14, 

-.13, 

-.12, 

-.10, 

-.09, 

-.10, 

-.09, 

-.07. 

-■05, 

-.04, 

+ 

-.03, 

.00, 

.01. 

-03, 

•04. 

•04, 

.04, 

.05, 

• 05, 

.05, 

+ 

• 04, 

•04, 

.04. 

•05. 

.05. 

.07. 

.08, 

.09, 

.10, 

.11, 

+ 

•10, 

.10, 

•09, 

.10, 

.10, 

•09, 

•09, 

.09, 

•09, 

.08, 

+ 

.07, 

•07, 

•07, 

•09, 

.09. 

.08. 

.08, 

.08, 

.08, 

.08, 

+ 

.06, 

-04, 

■oo. 

-.03, 

-.07, 

-.11. 

-.16, 

-.22, 

-.29, 

-.41/ 


BEGIN EXECUTION 


I F( ISTART.EQ.O) GO TO 30 

ISTART=0 

DO 10 I FREQ=1 ,4 

I F( I FREQ.EQ.2) GO TO 10 

XFAC=( 1 . -CHI ( 1 , I FREO)*CH 1(2,1 FREQ) )*( 1 . -DELTA< I FREQ) ) 
AW( I FREQ)=(1 .+CHI ( 1 , IFREQ))/XFAC 
AHV( I FREQ)=-CHI ( 1 , I FREQ)*(1 .+CHI (2, IFREQ) )/XFAC 
AOV(IFREO)=(1.-AW(IFREQ)-AHV(IFREQ»*2.7 
AHH( I FREQ)=( 1 .+CHI (2, I FREQ))/XFAC 
AVH( I FREQ)=-CHI (2, 1 FREQ)*(1 .+CHI (1 , 1 FREQ) )/XFAC 
AOH ( I FREQ)=( 1 .-AHH( 1 FREO)-AVH( IFREO)>*2.7 
10 CONTINUE 
C 

DO 20 ICEL-1,64 

SBIAS1 < I CEL )=1 . -SBIAS1 ( ICEl)/AVGTA(1 ) 
SBIAS2(ICEL)-1.-SBIAS2(ICEL)/AVGTA(2) 

SBIAS3C ICEL)=1 . -SBI AS3( ICEL)/AVGTA(3) 

SBIAS4( ICEL)=1 . -SBI AS4( ICEL )/AVGTA(4) 

SBI AS5( ICEL)=1 . -SBI AS5( ICEL)/AVGTA(5) 

20 CONTINUE 
C 

30 CONTINUE 
C 

DO 100 ICEL=1 ,64 

TA19V=TALO(1 , ICEL)*SBIAS1 ( ICEL)-ABI AS(1 ) 

TA19H=TALO(2, ICEL)*SBIAS2(ICEL)-ABIAS(2) 

TA22V=TALO(3, I CEL )*SB I AS3( I CEL )-ABI AS(3 ) 

TA37V=TALO(4, I CEL )*SB I AS4< ICEL )-ABI AS(4 ) 

TA37H=TALO( 5 , I CEL )*SB I AS5< ICEL )-ABI AS(5 ) 

C 

TALO( 1 , 1 CEL )=AW( 1 )*TA19V*AHV( 1 )*TA19H+A0V( 1 ) 

TAL0<2, IcEL)=AHH(1 )*TA1 9H*AVH ( 1 )*TA1 9v>A0H ( 1 ) 

TALO<3, ICEL)=1 .01993*TA22V+1 .994 

TALO(4, ICEL)=AW(3)*TA37V+AHV(3)*TA37H+AOV(3) 

TALO( 5 , 1 CE L )=AH H ( 3 )*TA37H+AVH (3 )*TA37V+AOH ( 3 ) 

100 CONTINUE 
C 

IFO85CHZ.EQ.0) RETURN 
C 

DO 200 ICEL=1 ,128 

TB85V=AW(4 )*ATAH I ( 1 , 1 CEL >+AHV(4 )*ATAH 1(2,1 CEL )+AOV(4) 
TB85H=AHH (4 )*ATAH 1(2,1 CEL )+AVH ( 4 )*ATAH I ( 1 , 1 CEL )+AOH (4 ) 
ATAHK1 , ICEL)=TB85V 
ATAHI(2, ICEL)=TB85H 

TB85V=AW(4)*BTAH 1(1,1 CEL )+AHV(4 )*BTAH 1(2,1 CEL )+AOV(4 ) 
TB85H=AHH(4)*BTAHI (2, ICEL)+AVH(4)*BTAHI(1 , ICEL)+AOH(4) 
BTAHK1 , ICEL)=TB85V 
BTAHI (2, ICEL)=TB85H 
200 CONTINUE 
RETURN 
END 
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SUBROUTINE FDTB00( I85GHZ) 


THIS SUBROUTINE CONVERTS ANTENNA TEMPS. TO BRIGHTNESS TEMPS. 
ALONG-SCAN AND OFFSET CORRECTIONS SET TO ZERO 

REAL*4 DELTA(4),CHI (2,4),AVGTA(5),ABIAS(5) 

REAL*4 SB1AS1(64),SBIAS2(64),SBIAS3<64),SBIAS4(64),SBIAS5(64) 
REAL*4 AW(4),AHV(4),A0V(4),AHH(4),AVH(4),A0H<4) 

SPECIFY COMMON /OUTDAT/ 

REAL*8 REV, XT I ME 

INTEGER*4 ITIME, ITIMSC, IVOLT, IAGC, ICOLDA, IHOTA, ICOLDB, IHOTB, IASCTM 
INTEGER*4 IATOIL, IBTOIL, I SAT 

REAL*4 XLATSC,XLONSC, ALTSC,THT, HLTEMP.RFTEMP, FRTEMP 
REAL*4 PERIOD, ASCLOC , ANG I NC , AX I S , ECC , ANGPER , SPACER 
REAL*4 ALAT , ALON , BLAT , BLON , TALO, ATAH I , BTAH I 
COMMON/OUTDAT/ REV, XT I ME , I T I ME , I T I MSC , XLATSC , XLONSC , ALTSC , THT , 

1 HLTEMP(3), IVOLT(2),RFTEMP, FRTEMP, IAGCC6), 

2 IASCTM, PERIOD, ASCLOC, ANGINC, AXIS, ECC, ANGPER, ISAT,SPACER(6), 

2 ICOLDA(5,7), IHOTA(5,7), ICOLDB(5,2), IHOTB<5,2), 

3 ALAT (128), ALON (128),BLAT(128), BLON (128), 

4 TALO( 5, 64), ATAH I (2, 128), BTAH I (2, 128), IAT0ILO28), IBT0ILO28) 

DATA INITIALIZATION 
DATA I START/ 1/ 

DATA DELTA/0.03199,0.02685,0.01434,0.01186/ 

DATA CHI/. 00379, .00525, .00983,0.0, .02136, .02664, .01387, .01967/ 

DATA AVGT A/ 190. 93, 130. 14,215.42,21 1.39,158.16/ 

DATA ABIAS/5*0./ 

DATA SBIAS1/64*0./ 

DATA SBIAS2/64*0./ 

DATA SBIAS3/64*0./ 

DATA SBIAS4/64*0./ 

DATA SBIAS5/64*0./ 

C 

C BEGIN EXECUTION 
C 

IF(ISTART.EO.O) GO TO 30 
ISTART=0 
DO 10 IFREQ=1,4 
IF(IFREQ.E0.2) GO TO 10 

XFAC=(1 .-CHI (1 , 1FREQ)*CHI (2, I FREO))*<1 .-DELTA(IFREQ)) 

AW( I FREO)=( 1 .+CHI (1 , 1 FREQ) )/XFAC 
AHV( I FREQ)=-CHI (1 , 1 FREQ)*( 1 .+CHI (2, 1 FREQ) )/XFAC 
AOV( I FREQ)=( 1 .-AW( I FREQ)-AHV( I FREQ))*2.7 
AHH( IFREQ)=(1 .+CHI (2, IFREQ))/XFAC 
AVH( I FREQ)=-CHI (2, 1 FREQ)*< 1 .+CHI ( 1 , 1 FREQ) )/XFAC 
AOH( I FREQ)=( 1 . — AHH C I FREQ)-AVH( I FREQ))*2.7 
10 CONTINUE 
C 

DO 20 ICEL=1 ,64 

SBIAS1 ( ICEL)=1 . -SBIAS1 ( ICEL )/AVGTA( 1 ) 
SBIAS2(ICEL)=1.-SBIAS2(ICEL)/AVGTA(2) 

SBIAS3( ICEL)=1 .-SBIAS3C ICEL )/AVGTA(3) 
SBIAS4(ICEL)=1.-SBIAS4(ICEL)/AVGTA(4) 

SB I AS5 ( I CEL )=1 . -SB I AS5 ( I CEL )/AVGTA<5 ) 

20 CONTINUE 
C 

30 CONTINUE 
C 
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DO 100 ICEL=1,64 

TA19V=TALO(1 , ICEL)*SBIAS1( ICEL)-ABIAS(1 ) 
TA19H=TALO(2,ICEL)*SBIAS2(ICEL)-ABIAS(2) 
TA22V=TALO(3,ICEL)*SBIAS3(ICEL)-ABIAS(3) 
TA37V=TALO(4,ICEL)*SBIAS4(ICEL)-ABIAS(4> 

TA37H=TALO(5, ICEL)*S8IAS5(1CEL)-ABIAS(5) 

C 

TALO( 1 , 1 CEL )=AW( 1 >*TA19V+AHV( 1 >*TA19H+A0V( 1 ) 

TALO(2, ICEL)=AHH(1 )*TA19H+AVH(1)*TA19V+A0H(1 ) 

TALO(3, 1 CEL >=1 .01993*TA22V+1 .994 
T ALO( 4 , 1 CE L ) =A W( 3 )*TA37V+AHV(3 )*TA37H+A0V(3 ) 

TAL0(5, ICEL)=AHH(3)*TA37H+AVH(3)*TA37V+A0H(3) 

100 CONTINUE 
C 

I F( I85GHZ.EQ.0) RETURN 
C 

DO 200 ICEL=1 , 128 

TB85V=AW(4)*ATAH 1(1,1 CEL )+AHV(4)*ATAH I (2 , 1 CEL )+AOV(4 ) 
TB85H=AHH(4)*ATAHI(2, ICEL)*AVH(4)*ATAHI(1 , ICEL)+AOH(4) 
ATAHKl , ICEL)=TB85V 
ATAHK2, ICEL)=TB85H 

TB85V=AW(4 )*BTAH 1(1,1 CEL)+AHV(4)*BTAH I (2 , 1 CEL )+AOV( 4) 
TB85H=AHH(4)*BTAHI(2, ICEL)+AVH(4)*BTAHI(1 , ICEL)+AOH(4) 
BTAHI (1 , ICEL)=TB85V 
BTAHI (2, ICEL)=TB85H 
200 CONTINUE 
RETURN 
END 
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